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CORRELATION, EVALUATION, AND EXTENSION OF LINEARIZED THEORIES FOR TIRE 

MOTION AND WHEEL SHIMMY 1 


By Robebt F. Smiley 

\ 


SUMMARY 

An evaluation is made of the. existing theories of linearized 
tire motion and wheel shimmy. It is demonstrated that most of 
the previously published theories represent varying degrees of 
approximation to a summary theory developed herein which is a 
minor modification of the basic theory of Von ScMippe and 
Dietrich. In most cases where strong differences exist between 
the previously published theories and the summary theory, the 
previously published theories are shown to possess certain 
deficiencies. 

A series of systematic approximations to the summary theory 
is developed for the treatment of problems too simple to merit the 
use of the complete summary theory, and procedures are dis- 
cussed for applying the summary theory and its systematic 
approximations to the shimmy of more complex landing-gear 
structures than have previously been considered. 

Comparisons of the existing experimental data with the 
predictions of the summary theory and the systematic approxi- 
mations provide a fair substantiation of the more detailed 
approximate theories. However, some discrepancies exist 
which may be due to tire hysteresis effects or other unknown 
influences. Thus, further work may be needed to explain these 
discrepancies. 

INTRODUCTION 

In the ground maneuvering of aircraft equipped with 
swiveling landing gears there sometimes arises the problem of 
violont oscillations or shimmy of the landing gear which may 
lead to failure of the gear. In the past this problem has been 
handled largely by means of various measures based on 
practical experience. However, this empirical approach has 
not proved entirely adequate. Moreover, for radically 
different types of complex flexible landing gears it is highly 
doubtful whether any empirical approach based purely on 
past experience could always safely and optimumly take into 
account all of the possible conditions which a landing gear v 
might bo subjected to in actual operation. 

A considerable amount of theoretical and experimental 
work on wheel shimmy has been done, mostly in the past .25 
years. (Most of the existing papers on this subject are 
listed in ref. 1, which also presents a historical discussion of 
the development of the wheel-shimmy problem.) However, 
most of these theoretical papers have not been correlated 
with each other or with the available experimental data, so 

* Supersedes NACA Technical Note 3632 by Robert F. Smiley, 1958, 


that essentially there exists at present a large number of at 
least superficially different theories of wheel shimmy and a 
fair amount of experimental data which has not been cor- 
related with many of these theories (refs. 2 to 23). 

The primary purpose of the present report is to clear up 
this partial confusion of theories by demonstrating that 
most of the previously published theories represent various 
approximations to one basic general linearized theory 
derived herein and that most of the previously published 
linearized theories which do not represent approximations 
to this general theory possess certain undesirable charac- 
teristics. This basic general theory, which is henceforth 
called the summary theory, is derived in such a manner that 
it makes use of and is compatible with the soundest features 
of practically all the previously published theories, insofar 
as this is possible at present; however, in the main this 
summary theory is a minor modification of the theory pro- 
posed by Yon Schlippe and Dietrich in references 3, 4, and 5. 

A second purpose of this report is to develop a series of 
systematic approximations to the summary theory suitable 
for use in the treatment of problems too simple to merit 
the use of the complete summary theory and to examine 
both these systematic approximations and the previously 
published theories to determine how these theories are 
related to the summary theory and how the predictions of 
these theories agree with the available experimental data. 

A final purpose of this report is to illustrate procedures 
for applying the summary theory and its approximations 
to complex types of flexible landing-gear structures. 

Although the primary purpose of this report is concerned 
with the wheel-shimmy problem, most of the material 
presented is directly applicable to the more general problem 
of the motion of elastic tires under arbitrary rolling condi- 
tions. Thus this material is pertinent to the study of 
veering-off or ground looping, ground handling, and cata- 
pulting stability of aircraft. 

The material in this report is arranged as follows. First, 
a detailed statement of the problem is given, together with 
a detailed outline of the manner in which it is treated herein. 
Then, after a brief discussion of the restrictions on the 
analysis, a linearized derivation is made for the general 
Von Schlippe-Dietrich type of kinematic equations govern- 
ing the motion of elastic tires rolling without skidding. 
This analysis proceeds essentially in accordance with the 
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theoretical analysis of Yon Schlippe and Dietrich except 
that the present analysis considers the subject of tire tilt 
in slightly greater detail. 

Next, the primary forces and moments acting on a rolling 
tire are discussed and used to establish the equations of 
motion for arbitrary rolling conditions. Then a systematic 
procedure is developed for forming approximations to the 
summary theory. 

The previously published theories are listed, discussed, 
and compared with the summary theory and these system- 
atic approximations. Finally, the application of the sum- 
mary and approximate theories to several simplified landing 
gears is discussed. The first example is chosen primarily 
to demonstrate the correlation between theory and experi- 
ment, the second example to demonstrate the correlation 
between the summary theory and its systematic approxima- 
tions, and the remaining example to illustrate the applica- 
tion of the theory to complex problems. 

Some of the material presented herein was submitted to 
the University of Virginia in partial fulfillment of the 
requirement for a Master of Aeronautical Engineering 
degree. 

SYMBOLS 


a 


a'—a—^Lh- 

■Ai,Ai,A3 

B u Bt 


C\ 

<h 

Cl 

<h 

di,d? 

D 

D( ) 

A( ) 


E§,E\, 

E 

f 

Eq,F t, . 

F w 

F„, 

F„ 

F v \ 


trail (perpendicular distance between ground- 
contact center point and swivel axis) 
sin k 
r 

coefficients defined by equations (115b) 
coefficients defined by equations (115b) 
lateral distance of center of pressure of vertical 
force from XZ-plane 

change in lateral distance of center of pressure 
of vertical force from ATZ-plane per unit of Xo 
change in lateral distance of center of pressure 
of vertical force from YZ-plane per unit of y 
distance from wheel center to center of gravity 
of swiveling parts of landing gear 
distance from center of gravity of swiveling 
parts of landing gear to swivel axis 
hysteresis constants used in equation (128c) 
tire parameter used by Bourcier de Carbon 
differential operator with respect to distance, 
d( ) 


dx 


or v 1 D t ( ) 


differential operator with respect to time, 

iror ) 

coefficients of linear differential equations 
energy dissipated per cycle 
frequency, v/2ir 

coefficients of linear differential equations 
lateral force due to hysteresis effects 
lateral inertia "force resulting from lateral de- 
formation of tire 

net lateral tire force acting on wheel 
net lateral structural force acting on wheel 
lateral force on tire due to lateral distortion of 
tire 


J-rr 

F, 

F, 

9 

h, 

I( ) 

/„ 


dfw 

I* 

j 

J 

k 

K 

K a 

AK a 

A K a{ 

A K a] 

K y 

Kk 

Affx 


A K X( 
A K h 


lo,li, • ■ • 

L 

m 

m, 

m u 

mi 

M 0 

M„ 

M* 

M* 

M zi 

■M„ 

M ta 


lateral force on tire due to lateral tilt of tiro 
vertical load on tire 

lateral force on swiveling parts of landing gear 
due to landing-gear strut 
linear damping constant (Damping moment== 
9 Dit) 

half-length of tire-ground contact area 


imaginary part of ( ) 

moment of inertia of the swiveling part of a 
landing gear about an axis parallel to the 
swivel axis and passing through the center of 
gravity of the swiveling part 
polar moment of inertia of wheel and tiro about 
an axis perpendicular to the wheel axle 
polar moment of inertia of tire (excluding solid 
wheel parts) 

total polar moment of inertia of wheel and tire 
moment of inertia of the swiveling part of a land- 
ing gear about the swivel axis 
excess of number of zeros over number of poles 
parameter in stability-determination plots 
(appendix C) 

lateral spring constant of landing-gear strut 
parameter used in appendix A 
torsional stiffness of tiro 

total effective change in the torsional stiffness 
due to the inertia effects 
effective change in torsional stiffness of tiro due 
to lateral acceleration of the 
change in the torsional stiffness duo to centrif- 
ugal forces 

lateral the force due to tilt per radian of tilt angle 
lateral stiffness of the 

total effective change in lateral stiffness of the 
due to the inertia effects 

effective change in lateral stiffness of tire duo to 
lateral acceleration of the 
change in lateral stiffness of the due to centrif- 
ugal forces 


the constants; 


, (nL+h)h H ~ i 
n\ 


relaxation length i 

mass of swiveling parts of landing gear 

mass of the 

mass of wheel including the 
mass of nonswiveling parts of landing gear 
constant friction-damping moment 
net structural- tilting moment acting on wheel 
center 

gyroscopic moment due to swiveling 
twisting moment due to hysteresis effects 
inertia moment resulting from lateral deforma- 
tion of the 

net structural swiveling moment acting on 
wheel center 

torsional moment on the due to twist of the 
gyroscopic moment due to tilting 
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gyroscopic moment due to lateral distortion of 
tire 

Mi torsional moment on the swiveling parts of the 

landing gear due to landing-gear strut and 
damper unit 

M^ damping moment about swivel axis 

N cornering power (lateral tire force per radian of 

yaw angle during steady yawed rolling for 
yaw angle approaching zero) 

72 1 parameter used in appendix A 

p complex roots of characteristic equations 

Pi,Pt functions defined in appendix D 

Pin,Pi'»,Pii ) Pi3,Pai,Pn functions defined in and after equa- 
tions (80) 

parameter used in appendix A 
r -free tire radius 

r, polar radius of gyration of tire 

r s vertical distance from wheel axle to ground 

radius of cross section of tire torus 
It tire parameter used l> 3 r Bourcier de Carbon 

R( ) real part of ( ) 

s circumferential coordinate on tire (fig. 1) 

S_ wave length, 2n/v l 

S tire parameter used by Bourcier de Carbon 

t time 

At time lag due to tire hysteresis 

T_ Moreland’s time-lag constant 

T tiro parameter used by Bourcier de Carbon 

T u Ti, . . . functions of D, correlating structural forces, 
moments, and deflections 

T a time-lag constant for hysteresis moment 

T\ time-lag constant for hysteresis force 

u,=C\F,(ffl t cos k— a)sin k 

u t i=c\F t (li cos k — a)sin k 

v rolling velocity 

w width of tire-ground contact area 

•Wi densit 3 r 

w,= (fiKy—aF t -\-CyF t sin K)sin k 

x horizontal distance parallel to mean direction of 

rolling motion 

X,T,Z space-fixed coordinate axes; the -X-axis is hori- 
zontal and parallel to the mean direction of 
rolling motion, the Z-axis is vertical, and the 
Taxis is perpendicular to the XZ-plane. 
The XT-plane is the ground plane, 
lateral distance of tire equator from XZ-plane 
vortical distance up from XT (ground) plane 
lateral deflection of tire equator from XZ-plane; 
subscript 0 refers to the center of the ground- 
contact area, 1 to the foremost point of the 
ground-contact area, 2 to the rearmost point 
of the ground-contact area, i to equator points 
, off the ground, and g to equator points on the 

ground 

a twist in tire, radians 

Aa half-width of twisting-moment — angular-deflec- 

tion hysteresis loop of tire 
7 lateral wheel tilt, radians 

lateral tire tilt resulting from lateral deforma- 
tion, radians 


V 

z 

VoPJuVi, 

VuVt 


VjVdjVijVzi 

Va>ViiVt 


y* 

Vi 

Va 

n\ 

e 

K 

A,Xo, Ai, 


AXo 


f 

p 

Pt= (ghK-y 


z l-{-^Lh 


Vi,cr 3 

T 



Subscripts: 


pneumatic caster, K a /N 

lateral deflection of center plane of wheel with 
respect to the XZ-plane; subscript 0 refers to 
the point corresponding to the center of the 
ground-contact area, 1 to the point cor- 
responding to the foremost point of the 
ground-contact area, 3 to the center point of 
the wheel, a to the point of attachment of the 
swiveling parts of the wheel to the swivel axis, 
i to wheel-plane points off the ground, and g 
to wheel-plane points on the ground 
inertia-force parameter (eq. (46)) 
inertia-moment parameter (eq. (49)) 
hysteresis-moment parameter (eqs. (62)) 
hysteresis-force parameter (eqs. (60)) 
angle of rotation of wheel about the vertical 
Z-axis, radians 

inclination of swivel axis, radians (fig. 5) 
lateral distortion of tire equator with respect to 
the solid parts of the wheel; subscript 0 refers 
to the center of the ground contact area, 1 to 
the foremost point of the ground contact 
area, 2 to the' rearmost point of the ground 
contact area, i to equator points off the 
ground, and g to equator points on the 
ground 

half-width of lateral-force — lateral-deflection 
hysteresis loop of tire 

circular frequenc 3 r of shimny motion, 2 irf or 
path frequency of shimmy' motion, vv~ l 
tiro tilt parameter (eq. (13)) 
spring constant for a linear restoring moment 
al<\ + acxFt T CyF t sin k) sin k 
tan k 
hr 

constants representing phase shift 
constant defined by equations (33) and (50) 

constant for gyroscopic moment ^ 


angle of rotation of wheel about the swivel 
axis, radians 

angular velocit 3 - of wheel about its axle ^ «;) 


max 


critical 

maximum 


STATEMENT OF THE PROBLEM AND GENERAL 
APPROACH 

The purpose of this section is to define specificalty tire 
problem considered in this report and to clarify further the 
correlation between the various parts of the report. 

The basic problem to be considered is the rolling motion 
and wheel s himm y of a rigid wheel equipped with an elastic 
tire, when the wheel is attached to some supporting structure 
such as a landing-gear strut. The motion of the rigid wheel 
can, of course, be completely described by six independent 
variables corresponding to the three degrees of freedom in 
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translation and rotation of the wheel. In addition to these 
six degrees of freedom, there exists a seventh degree of free- 
dom which is associated with the distortion of the elastic 
tire or the track of the tire on the ground which results from 
the application of a given motion to the rigid wheel. Thus, 
in general, the motion of a rigid wheel with an elastic tire 
represents a system of motion involving seven variables, and 
seven equations correlating these different variables are 
required to solve for the motion of a landing gear under 
arbitrary rolling conditions. Six of these equations will 
usually be the equations, expressing the sum of the forces or 
moments acting along each of the three principal coordinate 
axes; the seventh relation will be an equation, usually a 
kinematic equation, which correlates the tire distortion with 
the other variables. 

The present report is not concerned with all seven degrees 
of freedom. Most of the report is restricted to a considera- 
tion of cases of wheel motion in which the wheel is rolling 
at an approximately constant velocity v without braking, 
and consequently with constant angular velocity u, and 
where no strong vertical oscillations are involved. Thus, 
for example, effects of acceleration or deceleration, which 
are known to have at least some influence on the rolling 
motion (see, for example, the experimental evidence of ref. 
17)’ are not considered. Similarly, fore and aft oscillations 
of the wheel are excluded. 

When these three restrictions are applied, the seven- 
variable problem of a rolling wheel becomes reduced to the 
consideration of a system involving the following four degrees 
of freedom: (1) swiveling of the wheel about a vertical axis 
through the wheel center point, designated by the symbol 8; 
(2) lateral tilting of the wheel with respect to a vertical 
plane parallel to the direction of undisturbed motion, desig- 
nated by the symbol y; (3) lateral displacement of the wheel 
with respect to a space-fixed reference axis parallel to the 
direction of undisturbed motion, designated by the symbol y 
with various subscripts; and (4) lateral displacement of the 
tire footprint on the ground (which is a measure of the tire 
distortion), designated by the symbol y a . (These coordinates 
and their positive directions are illustrated in fig. 1.) 

In order to obtain four equations correlating these four 
variables 8, y, tj, and y a , the following procedure is used: 
After some remarks on general restrictions, a kinematic 
relation between the four variables is derived in the section 
entitled “Kinematic Relations for the Rolling Tire.” Next, 
the primary forces acting on the wheel from the ground, 
including wheel inertia forces, are discussed in the section 
entitled “Forces and Moments on the Wheel.” By utilizing 
these ground forces and moments, the four basic equations 
of motion for the wheel, including the kinematic equation, 
are set down in the section- entitled “Equations of Motion.” 

For many applications these equations of motion in their 
most general form are relatively complicated and, although 
they are by no means insolvable, it is profitable to simplify 
the equations for those problems which do not require the 
detailed equations of the summary theory. Therefore, a 
number of systematic approximations to the summary theory 
are formulated in the section entitled “Systematic Approxi- 
mations to the Summary Theory.” A second reason for 
establishing these systematic approximations lies in the fact 


that they furnish a framework for comparing the summary 
theory with the other existing theories of wheel motion, 
most of which are closely related to these systematic approx- 
imations. Such a comparison of the summary theory and 
its systematic approximations with the existing theories of 
wheel motion is carried out in the section entitled "Classifi- 
cation and Evaluation of Existing Theories.” 

In the last major section of this report the summary 
theory and its systematic approximations are applied to 
three illustrative types of landing-gear configurations which 
are chosen either to illustrate agreement between theory 
and experiment or to illustrate methods for applying the 
theory to complex problems of wheel shimmy. 

GENERAL RESTRICTIONS 

Before entering upon the detailed derivation of the equa- 
tions of motion, some further restrictions on the analysis 
should be discussed. First of all, the present report is 
limited almost exclusively to linearized theories. However, 
there is some question as to whether a linearized theory is 
sufficient to describe the important features of wheel shimmy. 
It appears at present that a linearized theory will provide at 
least a fair qualitative description of stability boundaries 
for shimmy and will indicate whether a given motion is 
stable or not. However, agreement between theory and 
experiment, presented in a subsequent section, is still not 
good enough quantitatively to warrant the conclusion that 
nonlinear effects can always be neglected or replaced by 
equivalent linear effects. 

Another limitation of. the linearized theory is that it doo9 
not permit calculation of the maximum steady-state shimmy 
amplitude for those steady-state self-excited shimmy motions 
which sometimes occur on actual landing gears. 

Although the preceding considerations suggest that non- 
linear effects in landing-gear motions possibly be of 
importance for, somp practical prob^ins, ( t^iq^p consideration 
is beyond the scope of the present report qnd henceforth 
only linearized theory is discussed. . The only concession to 
nonlinearity is made in appendix A, whjch presents a con- 
ventional approximate method for converting a nonlinear 
shimnly damper to an equivalent linear damper. It should, 
however, be noted that some attention has been given to the 
development of nonlinear tire-motion theory in references 
18 to 21. 

Another restriction arises in connection with the assump- 
tion adopted throughout this report that the finito width 
of the tire need not be taken into account in developing a 
tire-motion theory for single tires of conventional cross 
section. This assumption appears at present to bo at least 
partly justified on the basis of an experiment by Von Schlippe 
and Dietrich (ref. 3) ; on the other hand, since their investi- 
gation of this matter was extremely limited in scope, their 
experimental result may not be completely typical. Conse- 
quently, a more thorough evaluation of tire-width effects 
seems desirable. Some theoretical work, on this subject has 
been done by Von Schlippe and Dietrich (ref. 3) and later 
by Rotta (ref. 2), but the matter is beyond the scope of the 
present report. 

KINEMATIC RELATIONS FOR THE ROLLING TIRE 

In this section the kinematic equations for the motion of a 
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Fiodbb 1 . — Geometric relations for a rolling elastic tire. (For the sake of clarity parts (b) and (c) of this figure are drawn to different scales. 

The positive X-axis is the direction of undisturbed motion.) 


rolling tilted elastic tire without skidding are derived in 
accordance with the theoretical analysis of Yon Schlippe 
and Dietrich (ref. 3 or 4). This derivation differs only 
slightly from that analysis in that it omits some refinements 
of the theory which are necessary for very wide tires and it 
includes some influences of tilting of the the in more detail. 
While the modifications that are made in regard to tilt may 
not necessarily be of practical importance in most cases, 
they may bo of interest in a few problems. 

Specifically, the object of this section is to obtain a relation 
correlating the absolute lateral deflection of the center point 
of the the ground-contact area y 0 with the corresponding 
wheel coordinates of lateral deflection jj (for example, ij 0 or 
i? 3 ), swivel angle 6, and tilt y. (See fig. 1.) First, some 
geometric relations are set down and some background 
information regarding the distortion is discussed. Then this 
information is utilized to obtain a kinematic relation between 
the lateral deflection of the the center line or equator at the 
forward edge of the ground-contact area y t and the coordi- 
nates » 7 , y, and 6. Next, a kinematic relation between the 
lateral deflections of the the equator at the center and for- 
ward edge of the ground-contact area (designated y 0 and y u 


respectively) is established. These two relations are com- 
bined to obtain a basic kinematic equation correlating y 0 
with 77, 7, and 6. 

The derivation of these kinematic relations is based upon 
the following physical concept: As a tire moves forward, the 
tire material on the circumference just ahead of the ground- 
contact area is laid down or developed on the ground without 
skidding and becomes the new forward portion of the ground- 
contact area, so that the track of the tire is completely 
determined by the lateral-distortion coordinate of the fore- 
most ground-contact point y x and the slope of the distorted 
center line or equator of the the at that point. 

GEOMETRIC RELATIONS 

The primary geometric quantities involved in the problem 
of a rolling the are shown in figure 1, which gives an instan- 
taneous view of a distorted the with respect to an arbitrary 
space-fixed XYZ coordinate system, the X-axis being 
horizontal and parallel to the mean direction of wheel mot ion, 
the Z-axis being perpendicular to the ground, and the 
F-axis being perpendicular to the X- and Z-axes. Parts 
(a) and (b) of this figure represent side and bottom views, 
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respectively, of a rolling wheel that is swiveled and tilted. 
For the sake of clarity, part (c) of this figure, which shows 
an end view of the rolling tire, has been drawn to a different 
scale from part (b) and represents the unswiveled condition. 
In discussing the geometric quantities, the following termi- 
nology and symbols are used: The wheel center plane is the 
plane of symmetry of the wheel perpendicular to the wheel 
axle. The tire center line or equator comprises the tire 
points which on the undistorted tire are located at the 
intersection of the tire outer circumference with the wheel 
center plane; under the action of moments and lateral forces 
these tire points are deflected laterally by an amount X with 
respect to the wheel center plane. The symbol \ f designates 
the lateral deflection of tire equator points which are not in 
contact with the ground and \ t .designates the lateral 
deflection of points which are in contact with the ground. 
The point at the center of the ground-contact, area is desig- 
nated by Xo. 

The lateral distance of the wheel plane from an arbitrary 
space-fixed XZ-plane is designated by 7j, for points off the 
ground at a vertical height z and by for points on the 
ground. The lateral distances of tire-equator points from 
this XZ-plane are similarly designated by y { and y t . The 
difference between y and n is the lateral distortion X of the 


tire, or 




5 ? 

1 

£5 

II 

£ 

(1) 

and 




\ — Vt Vt 

(2) 


The tire contacts the ground in a finite area having a 
finite width and a length 2 h. The width of this area is 
assumed to be negligibly small; that is, the ground-contact 
area is assumed to be reduced to a ground-contact line. 
The foremost ground-contact point (in the direction of 
motion) is designated by the subscript 1, the rearmost 
point by the subscript 2, and the center point by the subscript 
0. Except for braking and accelerating effects, the center 
point 0 has approximately the same horizontal ^-coordinate 
as the wheel axle. 

Distances about the tire equator or circumference are 
measured in terms of the circumferential coordinate s whose 
origin is taken at the point 0. 

The wheel is assumed t& move at constant velocity v 
approximately in the direction of the X-axis. The wheel is 
laterally inclined with respect to the vertical Z-axis by the 
tilt angle 7 and is swiveled with respect to the XZ-plane by 
the swivel angle 0. Both tilt and swivel angles are assumed 
to be small; that is, cos 0=»cos 7 = 1, sin B~d, and sin 7=7. 

The center point of the wheel axle is located at a vertical 
distance r 3 from the XT (ground) plane, a lateral distance 
r 3 7 from the intersection of the wheel plane and the XT-plane, 
and a lateral distance % from the XZ-plane, where 

V 3 =vo—r 3 y ( 3 ) 

TIRE DISTORTION 

This section contains a short discussion of the features of 
tire distortion which are pertinent to the derivation of the 
basic kinematic relations of this report. 


Experimental and theoretical considerations (for example, 
see refs. 3 and 2 , respectively) indicate that, if the tiro 
equator in the ground-contact region is subjected to arbitrary 
lateral distortion, the lateral distortion of the tire equator 
off the ground X t tends to die out as an exponentially decaying 
function of the circumferential displacement s (for example, 
see fig. 2(a)). Thus, near tire point 1 off the ground the tiro 
distortion .will tend to approach the pattern described by 
the equation 

X,=V-'^ ( 4 ) 

and a similar equation will apply near tire point 2. The 
exponential constant I is a tire characteristic having the 
dimension of length and is called the relaxation length. 
The relaxation length near point 2 is not necessarily exactly 
the same as that near point 1; however, since the former 
relaxation length will not be used in this report in any 
critical calculations, this difference will not be taken into 
account. 

In regard to the accuracy of equation ( 4 ) very near point 
1, it should be emphasized that this exponential variation is 
only an expression of the equilibrium condition which the 
tire-equator distortion would reach in the absence of any 
restraints. However, it is obvious that conditions exist for 
which this distortion curve cannot' be completely exponential 
in form. For example, for the case of pure lateral deflection 
of a stationary tire, the tire equator in the ground-contact 
zone is (neglecting skidding) a straight line parallel to the 
wheel center plane and extending from point 1 to point 2 
(see solid lines in fig. 2 (b)). Consequently, tho existence 
of an exponential curve just to the right of point 1, and 
including point 1, would imply the existence of a sharp bond 
in the tire at point 1 such as is indicated in figure 2 (a). 
Since a sharp bend is impossible because of finite tire stiff- 
ness, it follows that, in general, on a stationary tire tho 
exponential variation given by equation ( 4 ) cannot bo valid 
close to point 1 . However, experimental evidence indicates 
that beyond a short transition region ahead of point 1 tho 
tire-equator distortion curve does have an essentially ex- 
ponential character (see solid lines in fig. 2 (b)). As tho 
wheel rolls ahead the nonexponential transition region of the 
tire equator is laid down or developed on the ground as it 
passes into the ground-contact zone, and the more nearly 
exponential part of the equator curve moves down toward 
the ground (see dashed lines in fig. 2 (b)) and is eventually 
developed on the ground, so that after rolling a short distanco 
from rest and during normal rolling conditions (fig. 2 (c)) tho 
tire-equator distortion at the front ond of tho tiro can 
approach the assumed exponential variation of equation ( 4 ). 

At the rear end of the tire the equator distortion curve 
during rolling does not so closely approximate an exponential 
variation, since at the rear end there is no process of laying 
down or development such as is responsible for the exponen- 
tial variation at the front end. However, since the rearward 
section of the tire equator is not used in any critical calcula- 
tions in this report, its equator curve is also, for simplicity, 
assumed to be exponential. 

If equation ( 4 ) is accepted as the basic equation for 
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(a) Assumed theoretical shape of tire equator distortion for a stationary 

twisted tire. t>=0. 

(b) Actual shape of tire equator distortion for an untwisted tire at rest 

(solid lines) and just after starting to roil (dashed lines). 

(o) Actual shape of tire equator distortion for a rolling tire. 

Figure 2. — Tire equator distortion. 

tire-equator lateral distortion near point 1 under rolling 
conditions, the total lateral displacement of the tire from 
the XZ-plano in this region can, by use of equation (1), be 
written in the form 

*-h 

yt=Vi+he L (5) 

Substituting the geometric relation y t = y,— yz (see fig. 1) 
into equation (5) gives 

t—h 

y i =y,—yz+'Kie L (6) 

KINEMATIC EQUATION 

By making use of the physical concepts discussed pre- 
viously, together with equation (6), it is now possible to 
establish as follows the basic differential equation relating 
the tire deflection at the center of the ground-contact area 
y 0 with the wheel coordinates y, 8, and y. 

There is assumed to be perfect adhesion between tire and 
ground, that is, no skidding. As the tire rolls forward 
(arbitrarily swiveling, tilting, and moving laterally) a dis- 
tance dx, a new element of the tire of circumferential length 
ds above and in front of point 1 is laid down or developed on 
the ground. This tire element, before being laid down on 
the ground, had the lateral-distortion variation given by 
equation (6). This equation, after differentiation with 


respect to s, yields for a given instantaneous position of the 
tire the following rate of change of distortion: 


$—h 


dy t dy t dz 1 a ( 7 ) 

ds~~ds 1 ds L Xie 

At point 1, where s=h and y t =y u 

< 8 > 

The term (jfe) is simply the sine of the angle between the 

ground and the tire equator at point 1. (See fig. 1.) Just 
to the left of point 1 the tire is flattened on the ground, or 

dz / dz\ 

^=0. If were not zero, the tire would have to have 

a sharp bend at point 1. However, because of the finite 
bending stiffness of an actual tire, a sharp bend is impossible; 

thus =0 and equation (8) reduces to 

(fH&K* » 

Further, since == ®’ s a h° r i zon t a l coordinate near 

point 1. The rate of change of wheel lateral displacement y t 
with respect to the horizontal coordinate x at any given 
instant is just the swivel angle 0; hence 


(dy A - 1, 

\dsJr L 1 


( 10 ) 


If the the is assumed to have no sharp bend at point 1, 

(jlsO = (^d^) P 0 ^- Then, since is the slope 

of the the equator on the ground at point 1 and since no 
skidding is assumed to exist, this slope must coincide with 

the track of the rolling the on the ground, which is 

Thus, ' 

dyi. 1 % 
fc~ d L ** 

or, if differentiation with respect to a: is designated by the 

d 

operator and the terms are rearranged, 


L Dyi=L0 — ( 11 ) 

(Alternate derivations of this equation are presented in refs. 3 
and 4.) A slightly more convenient form of equation (11) is 
obtained by substitution of the geometric relations Aj =y\—y\ 
and y l =y a - s rh8 (see fig. 1) to give 


(1 +£ D)yi=yo-\- ( L-\-h)8 


( 12 ) 
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Equation ( 12 ) is the basic equation for y x previously ob- 
tained by Von Schlippe and Dietrich (refs. 3 and 4). It 
should be noted, however, that no tilt terms appear in the 
equation. Although it is not known whether the effect of 
tilt on the validity of equation' ( 12 ) is important, in view of 
the present lack of a reliable method for taking this tilt 
effect into account, the following argument is presented to 
afford at least a crude approach to the problem. 

Equation ( 8 ) contains the tilt term y which was set 


equal to zero on the grounds that the factor 


is zero 


because of the finite bending stiffness of the tire. (See fig. 1 .) 
On the other hand, if it is assumed that the bending stiffness 
of the tire is zero and if radial tire distortion is neglected, 


(I),- 


be equal to h/r and the tilt term yhjr will enter into 


equations (9) to (12). As a somewhat questionable approxi- 
mation, it will now be assumed that a term of this type, but 
smaller by the reduction factor « 1 , should appear in the 
differential equations (9) to ( 12 ). Equation ( 12 ) then 
becomes 


(1 +L D) yi = Vo +'0(L+h )-& 


(13) 


A similar tilt term was derived by Greidanus (ref. 7) on the 
basis of a slightly different argument. (Greidanus’ term is 
discussed in a subsequent section of this report.) However, 
apparently no other detailed tire-motion theory has included 
such a term. 

Equation (13) is the fundamental kinematic relation for 
tire point 1 . The kinematic relations for points 0 and 2 are 
determined by the condition of perfect adhesion between tire 
and ground. During the rolling process each tire circumfer- 
ential element first contacts the ground at point 1 , later 
proceeds to point 0 and then to point 2 , after which it leaves 
the ground. Consequently, with perfect adhesion each tire 
element at point 2 has the same lateral deflection that it had 
when it entered the contact zone at point 1 a distance 2 h ago ; 
that is, 


y*{x)=y x (x—2h) 

Similarly the kinematic relation for point 0 is 
Vo(x)=yi (x—h) 


(14) 


or 


Vi(x)=y a (x+h) 


(15) 


Finally, by combining equations (13), (15), and (3) the 
equation 

[1 +L D}y,(x+h)=r, Q (x)+[L+h]e(x)~^(x) 

=, 3 W+[i+A]e(x)+[r 3 -^]y(x) (16) 


is obtained. This is the basic kinematic equation correlating 
the tire lateral deflection y 0 under rolling conditions with tho 
swivel angle 0 , the wheel lateral displacement % or %, and 
the lateral tilt y for arbitrarily applied variations of 0, ij 0 or 
7 j 3 , and y. However, this transcendental form of the kine- 
matic equation is not the most convenient form for some 
purposes in this report. In particular it is expedient to 
remove the transcendental expression from equation (16) 
by use of a series expansion. 

SERIES EXPANSION OF KINEMATIC EQUATION 

The expression y<>(x+h), after expansion in a Taylor series, 
gives 

2/o (x+h)=y a (x)+h Dy 0 (x)+y Lr-y 0 (*)+ . . . 


= 2 /o (x )+2 £>" 2/0 (x) 


(17a) 


where the operator D n represents An alternative form 
of equation (17a) which is useful later is 


y 0 (x+h)=e’ ,D 2/0 (x) 


(17b) 


since the infinite series in hD is tho series expansion of the 
exponential function. A third useful foim of this equation 
is obtained by expressing equation (17b) in terms of a timo 
derivative D t instead of the spaco derivative D. Since it 
is assumed throughout this report that the rolling velocity 
v is constant, tho correlation between these two derivatives 
is given by the equation 

■n / \ d{ ) dx d(_ ) d( ) T-,/ 

* 

and hence equation (17b) can also be written in the form 

yo(x+h)=e bf ~' 1), yo(x) (17c) 

Differentiation of equations (17) gives the result 


Dy 0 (x+h)=Dy 0 {x)+h D x y 0 (x) + 2 h t D 3 y 0 {x)+ . . . 


&•*<*> 


(18a) 


or 


Dy 0 {x+h)=Dd ,t> yo{x)^v- l D t d ,,r ‘ D 'y^x) (18b) 

Substitution of equations (17) and (18) into equation (16) 
gives after rearrangement, with y 0 (x) written simply as y n 
and similar treatment of 77 0 , 173 , 0 , and y: 

ytr\-lid ~~ l=Vz J rhS J r(xz~ 


= (l+fi D-\-li Z) 2 -}- . . .)y 0 
=(l+iblnD n '^i/o (19a) 
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where 

U={2L+h)\ 
l.= (nL+h) ^ 
or 

vo+lid—^~~ t =, ?3+ M+(W— 

= (1 +L D)e XD yo 

=(l+Lv- 1 D l )e hc ~' D ‘y 0 (19b) 

Equations (19a) and (19b) are alternative forms of the basic 
kinematic equation (16) which axe useful in subsequent 
sections of this report. 

This concludes the derivation and discussion of the basic 
kinematic equation correlating the lateral tire deflection 
y 0 -with the wheel coordinates 6, y 0 , and y. Next, attention 
will bo directed to the relationships existing between these 
coordinates and the forces and moments acting on the wheel. 

FORCES AND MOMENTS ON THE WHEEL 

In this section the primary forces and moments acting on 
a rolling wheel are discussed and, where possible, equations 
are sot down for these quantities. These equations are then 
utilized in later sections, together with the preceding kine- 
matic equation, to establish the equations of motion for a 
rolling wheel. 

The forces and moments considered fall into five general 
categories: elastic forces and moments due to tire distortion, 
gyroscopic moments, tire inertia forces and moments, 
hysteresis forces and moments, and structural forces and 
momonts. 

Throughout this discussion, forces along the coordinate 
axes are considered positive if they tend to move the wheel 
in the positive directions of the coordinate axes; moments 
about the coordinate axes X, 7, and Z or other parallel 
axes are considered positive if they tend to produce wheel 
rotation from the positive F-axis toward the positive Z-axis, 
from the positive Z-axis toward the positive X-axis, and 
from the positive X-axis toward the positive F-axis, respec- 
tively. 

ELASTIC FORCES AND MOMENTS DUE TO TIRE 
DISTORTION 

LATERAL ELASTIC FORCE 

The lateral elasticity properties of a tire will be considered 
first. If a static untilted tire is laterally deflected at its 
base with respect to its rim by a lateral force F v \, it produces 
an equal spring reaction force roughly proportional to the 
mean lateral distortion 'k meaH) or, inversely, a lateral tire dis- 
tortion X mean creates a proportional ground force F v \. If the 
lateral distortion of the center of the ground-contact line Xo 
is taken as the mean distortion, then the elastic ground force 
is 


where K\ is the lateral spring constant or side stiffness of the 
tire. This relation is used by most investigators. However, 
in references 2 to 5 a slightly different expression is used. 

In these references the mean lateral distortion of the tire is 
defined as the average of the distortions at the leading-edge 
and trailing-edge points of the ground-contact area (points 
1 and 2). The resulting equation for F v \ is 

X r x= ? Xx(X 1 +X 2 ) (21) 

instead of equation (20). The true equation for F„\ is prob- 
ably more complicated than either of these two equations; 
however, since no plausible means of obtaining a better 
equation is available, it appears advisable to select one of the 
above equations for use in this report. Equation (21) may 
be slightly the better equation for a few special cases of wheel 
motion, but equation (20) is much simpler to work with, and 
in most cases of wheel motion it makes little difference 
which of the two equations is used. Therefore, for the sake 
of simplicity equation (20) is adopted hereinafter as the 
basic equation for the lateral force on a wheel due to lateral 
deformation of the tire. 

TORSIONAL ELASTIC MOMENT 

The torsional elasticity properties of a tire will be con- 
sidered next. If a tire is twisted on the ground about a 
vertical axis through an anglo a, there arises a restoring 
ground moment that is roughly linearly proportional to 
the twist: 

M, a =K«a (22) 

The tire twist a is equal to the mean angle between the track 
of the tire on the ground and the wheel plane; that is, 
a =Dy m>an —e. Taking the value of Dy wan as Dy 0 gives 

a=I>y 0 -e (23) 

and thus 

M ia =K a ( Dy 0 — 6 ) =K a (y~ 1 D,y 0 —0) (24) 

Most investigators of tire motion use this relation. However, 
in references 2 to 5 the mean angle is taken equal to 
(Xj— X 2 )/2 /i and thus the moment equation 

(Xi— X 2 ) ' (25) 

is obtained, which leads to relatively more complicated 

equations of motion than does equation (24). Since there 

is no strong reason for believing equation (25) to be a 

significant improvement over the simpler equation (24), the 
latter is used in the analysis of this report. 

Melzer (ref. 10) has used the less accurate relation that the 
moment due to tire twist is 

M ta =—Kj) (26) 

which implies the relation ff^>Dy 0 . (See eq. (24).) Since 

this relation is not true in all practical cases, Melzer’s 
theory should be viewed with some caution. 

TILT ELASTIC FORCE 

If a tire is tilted from the vertical Z-axis by an angle V ' 


F t \—K\\ 1 =K\(yo— i?o)=Xx(y 0 — vs~ m) 


( 20 ) 
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without lateral distortion of the equator (A<>=0.), there arises 
a restoring ground lateral force that is approximately linearly 
proportional to the tilt -angle (e. g., see ref. 2); 

F„=-K y y (27) 

where K y is the constant of proportionality. Most authors 
(excepting Rotta in ref. 2) have not considered the effects of 
this force term although they have considered other effects 
of the same order of magnitude. 

VERTICAL-LOAD CENTER OF PRESSURE 

Under some circumstances the vertical load F z influences 
the wheel motion. In order to consider this influence it is 
necessary to know the location of the center of pressure of 
this force. In the ATZ-plane (fig. 1) this center of pressure 
lies approximately below the wheel axle in line with the point 
0. In the FZ-plane the center of pressure is shifted laterally 
from the intersection of wheel plane and ground i?o as a result 
of lateral distortion Xo and tilt y. As a first approximation, 
this shift may be taken as linearly dependent on Xo and y so 
that the lateral distance c of the center of pressure from the 
A r Z-plane becomes 

C= J?o+CxXo— -Cyy 
=cxy 0 + (1— Cx)tio— c y r 


2 



Figure 3. — Effective tire tilt due to lateral distortion of tiro. 


where I 1t is the moment of inertia of the tire (excluding tlio 
solid rim and axle) about the wheel axle. By using the 
relation D,( )=v D( ), equation (29) can also be expressed 
in the form 

M*=-I„v>^Dy x (30) 


=cxyo+ (1 — Cx)i73+ [(1 — Cx>3— Cy]y (28) 

where c\ and c y are constants. (The signs of the terms are 
chosen so that cx and c y are positive numbers.) 


GYROSCOPIC MOMENTS 


GYROSCOPIC MOMENT DUE TO LATERAL DISTORTION OF TIRE 


The origin of gyroscopic moments on a rolling untilted 
wheel with lateral distortion of the tire at the ground (fig. 3) 
is considered next. While the solid rim and axle parts of 
the wheel are untilted, lateral deformation of the elastic tire 
causes the tire, on the average, to be tilted with respect to 


the wheel center plane by an amount -yx= 


\>n 

r+r 3 ’ 


where r is 


the tire radius and ti is a correction factor which indicates 
the effective fraction of the total tire mass that is tilted at 
this angle. Kantrowitz (ref. 8), apparently the only in- 
vestigator who has considered this at least theoretically 
significant factor, has suggested that t, = 1/2. This tilting 

action produces an angular velocity D,yx=^~^-j where D, 

' f-|-r3 

indicates differentiation with respect to time. This angular 
velocity, together with the rotational velocity of the tire &>, 
produces a gyroscopic moment about the Z-axis of magnitude 


where the ratio r/co is, to a good enough approximation for 
this secondary term, equal to the tire radius r. Then, 
substituting for y\ and u/v in equation (30) gives 

< 31) 

For later convenience, the result can be expressed in several 
alternate abbreviated forms: 


M t \= — rv t D\o= — TFD(y 0 — Vo ) =—tv D,(y 0 — ns— m) (32) 


where 


_ zila 

r(r+r 3 ) 


(33) 


Another method for deriving an expression for r is discussed 
in a,subs£quent. section. 

GYROSCOPIC MOMENT DUE TO TILTING OF WHEEL 

If the entire wheel structure tilts at an angular volocity 
D,y, another gyroscopic moment arises of magnitude 

M rr =—I ru fo Fa « -Zf-D" (34) 


in addition to the term of equation (29). Here /„„ is the 
total polar moment of inertia of the wheel (including the 
tire) about its axle. 


M zk =—I tt w D t y\ 


(29) 
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GYROSCOPIC MOMENT DUE TO SWIVELING OF WHEEL 

If the wheel swivels at an angular velocity Dfi, a tilting 
gyroscopic moment also arises of magnitude 


M#= -I suM Dfi « Dfi (35) 

T 

TIRE INERTIA FORCES AND MOMENTS 

This section is concerned with an examination of the in- 
fluence of tire inertia forces and moments on a wheel rolling 
at high speeds. Two types of such inertia effects are evalu- 
ated now in separate sections: inertia forces and moments 
associated with lateral distortion and twisting of the tire, and 
centrifugal forces and moments. Then the overall effects 
of these two types of inertia forces and moments are con- ' 
sidered in another section. 

INERTIA FORCES AND MOMENTS DUE TO LATERAL DISTORTION OF TIRE 

At high rolling or shimmy velocities, tire inertia forces and 
moments arise which are proportional to the relative ac- 
celerations of the different parts of the tire (including the" 
previously discussed gyroscopic moment due to tire lateral 
distortion, discussed here from a slightly different point of 
view). A rough estimate of these forces and moments can 
bo made as follows: One-third of the total mass of the tire 
m t is assumed to be located on the periphery of the tire and 
to be subjected to the same accelerations with respect to the 
, wheel hub as are tire particles on the equator line, while the 
remaining tiro mass is assumed to be substantially undis- 
turbed. The “active” mass of the tire per unit circumfer- 
ential length is then mjGirr. The lateral acceleration of tire 
particles on the right-hand side of the tire and off the ground 
in figure 1 (a) will be considered first. The lateral distortion 
of the tire in this region is given by equation (4). The 
lateral relative velocity of a tire particle,, obtained by differ- 
entiating this quantity with respect to time, is 

D t \ i =^D t \ i — -j-D t s^e L 

The quantity Ds, which represents the peripheral velocity 
of tire particles with respect to the wheel axle, is approxi- 
mately equal to the negative of the rolling velocity v, so that 
the velocity expression becomes 

Differentiation of this result to give the relative acceleration 
of the tire particles yields the result 

A%=(#A i+|° Z^-f^x,)*" ^ 

The corresponding inertia force A F for this part of the tire 


D[\i — ^D t \i 


is obtained by integrating the product of this acceleration 
and the active mass per unit length to obtain the force 

term — I D?\ t ds. Evaluation of this integral, 

after replacing the upper limit by infinity for simplification 
of the result (which introduces no significant error because 
of the rapidly decaying exponential function in DfiXi), yields 

A^= AX, +^X : ) (36) 


The corresponding inertia moment AM is given by the 
expression 


m, r*~ rr 
6rr J »-h 


r sin <p Dfi\ t ds 


where r sin <p is the moment arm (see fig. 1 (a)), « is related 
to <p by the relation s— h=r(<p— <pi), and y>i=sin -1 —■ There- 
fore, the moment integral may be written in terms of p in 
the form 

AM== — £ P r Sbl p(A*Xi+f r d v 


Evaluation of this integral, after replacing the upper limit 
by infinity (which introduces no appreciable error), yields 
the expression 


m t rL 




m ‘ M&+F) ' (WX.+fu,X,+Yx,) (37) 


In a similar manner, for tire particles off the ground on 
the left-hand side of the tire in figure 1(a), the following 
expressions are obtained for the inertia force and moment: 

AF=~~^(^D, t X 3 — -j-D t \ 2 - (38) 


m,rL 
A M= 


6t KA 2 +ri) 



In these two expressions it is assumed, for reasons previously 
discussed, that the relaxation length L is the same for both 
sides of the tire. 

In obtaining the inertia forces and moments for tire par- 
ticles in the ground-contact area, it is recognized that in 
practically all cases where inertia forces are' important the 
ground-contact line is almost a straight line, so that the 
lateral distortion for tire particles in this region can be 
expressed fairly well by the equation 


X*=Xo+sa 
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The corresponding velocity and acceleration are 
D t \ t =Dt\i-{-8 D t ct — va 
D l i \ g =D i '\( > -\-s D, 2 a — 2d D,a 
The total inertia force for this region is then 


the dynamic case the relation between ground force and 
lateral distortion of the tire may be modified by the inertia 
effect. As a first approximation, it will be assumed that the 
modification of the ground force is proportional to tho inortia 
force, or 

F,\=K\\) — (4G) 


=-j^{D l ^-2vD l a) 
and the inertia moment is 


AM=-^ P"* s D?\ ds 

6TTTj t r,- h * . 


m,h 3 

9 7rr 


D, 2 * 


(40) 


(41) 


The total inertia force Fp,, obtained by summing the force 
■terms in equations (36), (38), and (40), can be stated 
conveniently in terms of Xo and a by using the relations 
7k]+Xj=2Xo and \ — \ 2 =2ha, which are valid for a sub- 
stantially straight ground-contact line. The result is 


Fn (li Xo ^ (42) 

Avliere li—L+h. Similarly, the total inertia moment is 



The significance of these inertia expressions will now be 
partly evaluated by considering the inertia force for sinu- 
soidal oscillations; that is, where X 0 =X 0>tai . sin vt and there- 
fore ZVX 0 =—r 2 Xo, so that equation (42) may be restated as 


F «—£r(ir 1 ^ <«> 


(An equation similar to eq. (44) has been derived by Mar- 
strand in ref. 20. Marstrand’s equation, however, is based 
on a cruder representation of the shape of the lateral dis- 
tortion of the tire.) 

In order to interpret the significance of the inertia force 
it is noted that the tire force quantity which is of importance 
for the subsequent analysis is the net tire force F^ acting 
on the wheel, which is equal to the sum of the ground force 
F v \ and the inertia force F ti : 

F„=F B x+F tt (45) 

Tor a static tire, F v \ was set equal to (eq. (20) ). In 


where ^ is a number whose absolute value will be loss than 
unity if the modification of the ground force duo to tho 
inertia force is less than the inertia force itself. 

After combining equations (44), (45), and (46), the follow- 
ing equation for the net tire force F tlt is obtained : 


Fv*— ^ttX+T )rt = ^2Cx— (1— ' >?y) ^ X° (47) 

From the form of equation (47) it can be seen that, insofar 
as the ratio of net tire, force to lateral deformation is con- 
cerned, the effect of the inertia force can be considered equiva- 
lent to a change in tire lateral stiffness AK\ { equal to 

AK h =-( l-^)g(x~^) (48) 


Similarly, from an examination of tho terms containing a 
in the inertia moment equation (43), it can be concluded 
that part of the effect of this inertia moment is to change tho 
tire torsional stiffness by an amount AK av winch iff defined 
by 


AK ai 


(1 

3tt 


-rh 

h+L-yJ 



hV 

_ 

L 2 - f-r 3 

\L L ^J 

“irj 


(49) 


where 77, is a number representing the . torsional stiffness 
similar to for the lateral stiffness. The remaining inertia- 
moment term in equation (43), which is proportional to 
£>,Xo, is simply the previously discusssed gyroscopic moment 
due to lateral tire distortion. By comparing this term "with 
equation (32) it is seen that the coefficient v may be expressed 
by the equation 


2m, r 
t— 


( h + L ^f-^) 

ML 3 +r*) 


(60) 


Equation (50) gives approximately the same result as equa- 
tion (33) with Kantrowitz’ assumption that r, =<i. Tho 

z 

discussion of the velocity range in which these stiffness 
changes are important is postponed until after the effects of 
centrifugal forces have been considered. 

EFFECTS OF CENTRIFUGAL FORCES 

Another inertia effect that may become significant at high 
speeds is produced by the centrifugal forces acting on the 
individual mass elements of the the. These centrifugal 
forces appear to increase the the stiffness, as ■will bo demon- 
strated by a crude analysis which gives a qualitative idea of 
this effect but which should not be regarded as possessing any 
strong quantitative merit. 

For the purpose of this estimate, one-half the mass of the 
the is assumed to be concentrated in the side walls and the 
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other half is assumed to be concentrated on the periphery. 

If the tire lateral and torsional stiffnesses K\ and K a are" 
assumed to be directly proportional to the tension in the side 
walls of the tire, there will he two sources of tire stiffness : 
inflation pressure, which produces a side-wall tension approx- 
imately equal to wp per unit circumferential distance (where 
w is the tire width), and centrifugal force, which produces the 


side-wall tension (sTt) (jr) corresponding to the pe- 


ripheral tire mass ~ m,. Thus the lateral stiffness of the 
tire may he expressed in the form 


Kxccbrt^wp-^rritV 1 

or, equivalently, as 




\taiic 


iirrhvp 


It is evident from this equation that centrifugal force increases 
the tire lateral stiffness by an amount A K\p 


A K\ 


^rhop 


(51) 


and the torsional stiffness by an amount A K aj : 


AK aJ 


irhtfK, 


a static 


i-trrhjop 


(52) 


SIGNIFICANCE OF TIRE INERTIA EFFECTS WITH RESPECT TO TIRE 
STIFFNESS 


The significance of the two just discussed tire inertia effects 
on the tire stiffness 'will now be considered. 

For the lateral stiffness of a tire, the effective change 
A K\ from the static value of K\ is obtained by adding the 
two increments given by equations (48) and (51). The 
resulting effective overall change in tire lateral stiffness as a 
function of rolling speed and shimmy frequency is 


A ^ (1 —V^mthv 2 (1 — 7 h )m^ 

s 37r r 3nr L i-rrhiyp 


(53) 


The fust term, which involves the shimmy frequency, 
appears to be small enough in comparison with K\ so that 
it can probably be neglected for most practical conditions. 
The last two terms have opposite signs if i?„<T and thus 
may represent two partly counterbalancing effects. The 
second term arises from the previous considerations of the 
lateral acceleration of tire particles and tends to reduce the 
effective lateral stiffness of the tire with increasing rolling 
velocity if 1. The last term arises from the previous 
considerations of centrifugal forces and tends to increase the 
lateral stiffness. These last two terms indicate that at high 
rolling speeds, if ,the tire stiffness may either drastically 
decrease or drastically increase, depending on which of the 
two terms is larger. However, both terms happen to be of 
the same order of magnitude and the derivations of both 


terms are based on concepts too crude to justify conclusions 
regarding which term is larger. Thus, the only conclusion 
that can be drawn is that at sufficiently high rolling speeds 
drastic changes in tire lateral stiffness may occur. Whether 
the stiffness increases or decreases can probably be settled 
only by experiment. 

In order to give some quantitative measure of the velocity 
at which these inertia effects become significant, some calcu- 
lations were made' to determine the velocity at which the 
magnitude of the second term in equation (53) becomes equal 
to /<x. By making use of the static tire data in reference 24 
for several modern aircraft tires and assuming that y u =0, 
it was found that this velocity averaged approximately 
400-\/r fps « 270 Vr mph, where r is expressed in feet. Similar 
estimates for the velocity at which the third term in equation 
(53) becomes equal to K\ yielded approximately this same 
velocity. Moreover, since this velocity is rather high com- 
pared with normal present-day landing speeds, the inertia 
effects on tire lateral stiffness considered here can probably 
usually be neglected. 

For the torsional stiffness of a tire, the overall effective 
change in torsional stiffness A K a due to tire inertia and 
centrifugal forces is obtained by adding the two increments 
given by equations (49) and (52). The result, is 


A K a 


(1—7 u)m t 

rh 
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^rhpp 


(54) 


This equation is parallel to equation (53) for the lateral 
stiffness, so that statements made previously concerning the 
lateral stiffness apply here also. 


OTHER INERTIA EFFECTS 

The preceding discussion suggests that the effects of tire 
inertia are to change tire stiffness at high speeds and to 
introduce a gyroscopic moment. However, it should be 
recognized that other inertia effects will come into play, 
probably at velocities close to those at which the previously 
mentioned inertia effects arise. For example, the basic 
kinematic equation depends on the assumption of an expo- 
nentially distorted tire-equator line corresponding to a 
definite “static” relaxation length. This assumption is valid 
(if it is valid at all) only when the elastic forces in the tire 
predominate over the inertia forces. Where inertia forces 
are strong in comparison with elastic forces, it is at least 
doubtful whether the relaxation length remains constant. 

Although there are undoubtedly other effects of tire inertia 
in addition to the ones discussed here, it appears probable 
that the importance of many tire inertia affects can be 
assessed by means of the following s umm ary statement: 
The major effects of tire inertia on the rolling motion appear 
to come into play at a velocity of an order of magnitude of 
400-yr fpss»270V?* mph where r is expressed in feet. For 
considerably smaller velocities, most inertia effects can 
probably be safely neglected; for velocities of this order of 


400104—0! 
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magnitude or higher, many of the basic assumptions of this 
report, and of most other papers on this subject, may be 
subject to considerable error. 

HYSTERESIS FORCES AND MOMENTS 

In addition to the 'forces and moments previously dis- 
cussed, certain damping forces and moments .arise as a 
consequence of the sometimes considerable hysteresis losses 
which arise in the distortion of elastic tires. Apparently 
the only significant attempt to deal with this hysteresis 
problem is reported by Von Schlippe and Dietrich in refer- 
ence 5. s This reference provides some valuable insight into 
the fundamental mechanism of the hysteresis process and 
presents an equation for the hysteresis moment acting, about 
the swivel axis of a shimmying wheel. However, even 
though use of this hysteresis-moment equation leads to good 
agreement between theoretical and experimental stability 
boundaries for a limited amount of experimental data (as is 
shown subsequently in the present report), some parts of 
the analysis seem so unrealistic that it is questionable whether 
much confidence can be placed in the final results of refer- 
ence 5. Apparently, the only other significant contribution 
to the hysteresis problem is provided by the analysis of 
Moreland in references 11 and 12. In these references tire 
hysteresis forces as such are not considered, but -the idea is 
introduced that a tire possesses a characteristic time-lag 
constant. In a subsequent section of the present report it 
is shown that this time-lag constant may be, at least in part, 
a consequence of hysteresis effects. However, the interpre- 
tation of Moreland’s time-lag constant as a hysteresis effect 
presents some questionable features that are also discussed 
subsequently. 

No completely satisfactory solution of the hysteresis prob- 
lem has been found yet. However, the following crude 
analysis of this problem offers another point of view with a 
few qualitative merits not possessed by the two previous 
analyses. 

Consider the case in which a standing tire is subjected to 
a periodical lateral deformation Xo of the form 


tire deformation \ lags behind the lateral ground force F v \ 
by a time increment At, where At is approximately equal to 
the ratio of AXq to the maximum slope of the curve of Xo 
plotted against time (which is 'K 0mnx v for the assumed varia- 
tion of Xo): 

At~^~ (66) 

As a first approximation for quantitatively including this 
time-lag concept in the present analysis, a tire is assumed to 
behave somewhat like a combination of a linear spring and 
a damper unit such as is indicated in sketch 3, where the 



spring constant corresponds to the previously discussed the 
lateral stiffness K\ and Cj is the coefficient of an equivalent 
linear damper. Inertia forces are neglected for the present 
argument. The differential equation for this system is 

jKxXo"h ®i D l \ 0 =F v \ (66) 

and its solution for the case of X 0 =X 0 „ WI sin vt gives a hystere- 
sis loop of the form indicated in sketch 2 where the time lag 
At becomes 


Under these conditions the lateral ground force F v \ on the 
tire is experimentally observed to vary with time in the 
manner indicated in sketch 1 and the corresponding varia- 
tion of lateral ground force F,\ with lateral tire distortion 
Xo, shown in sketch 2, appears in the form of a typical 
hysteresis loop. As can be seen from sketch 1, the lateral 






5 Although Von Schlippe and Dietrich considered hysteresis effects in an earlier papa 
(ret 4), this earlier analysis leads to some conclusions which are not in agreement with the 
results of the lata, more detailed analysis of reference 5. 


At 


tan“ l 


(67) 


After equating equations (65) and (57), can bo expressed 
isTx tan(AXo/Xo ) 

by the relation a 1 = — ■ — ■ so that equation (66) can 

then be written in the form 

F t \=K\Xf ) J t I'Va= AT x(XoT T\ D,Xo) (68) 


where F v „, the lateral force resulting from hysteresis effects, 
is 


F, h =K\T\ D,\o=K\J\v D \ o 

and where 

Ta — tan (AX°/X<w)_^ ' 

V V 

, AXo 
7ft=tanr — — 

"^atax 


(69) 


( 60 ) 
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With the same type of reasoning the hysteresis twisting 
moment is given by the equation 

M, b =K a T a D t a=K a T a v Da (61) 

(compare with eq. (59)), where 

T _ tan (Aa/Q _V« 

“ v v 

V (62) 

. Aa 
rj^tan 

It is now seen that the determination of the hysteresis 
force and moment from equations (59) and (61), respectively, 
depends on the determination of two quantities T\ and T a 
which have the dimension of time and which will be called 
time-lag constants. The quantity T\, in particular, can be 
considered closely analogous to Moreland’s time-lag constant. 
In the present case equations have been derived for the time- 
lag constants as functions of the three variables ip, ij a , and v. 
(See eqs. (60) and (62).) However, it should be clearly 
recognized that equations (60) and (62) are based in part on 
arguments valid only for a standing tire. These arguments 
may no longer be valid for a rolling tire, and oven if the idea 
of a time-lag constant is still valid, it is likely that the time- 
lag constants will not be adequately predicted by equations 
(60) and (62), particularly if the quantities Vv and rj, are 
evaluated from static hysteresis loops. Moreover, it is basi- 
cally unsound to assume that the hysteresis force is dependent 
only on the tire lateral distortion Xo and is independent of the 
tire twist a; actually the hysteresis force (and the hysteresis 
moment) will in general depend in a complex manner on both 
Xo and a, and even for a first approximation the interaction of 
these two variables cannot necessarily be neglected. Thus it 
appears that the preceding hysteresis equations are based on 
rather speculative and perhaps unsound assumptions, at 
least from a quantitative point of view, and for this reason 
these equations will not be incorporated into most of the 
derivations in subsequent parts of this report. On the other 
hand, the preceding derivation may be sufficiently plausible 
to give some idea of tho order of magnitude of hysteresis 
effects, particularly since Moreland has indicated in reference 
12 that his experimental data (mostly unpublished) demon- 
strates the existence of a time-lag effect in tire motion; con- 
sequently, in a few parts of this report some mention will be 
made of the consequences of introducing the hysteresis force 
and moment terms that have just been derived into a wheel- 
sliimmy analysis. 

STRUCTURAL FORCES AND MOMENTS 

The preceding discussion covers the major ground forces 
and moments and the gyroscopic moments acting on the 
wheel. In addition, forces and moments are exerted on the 
wheel by the supporting structure. These will be designated 
as F b , for the net structural force parallel to the T-axis, M x , 
for the 'net structural lateral tilting moment, and M„ for the 
not structural swiveling moment. These forces and moments 
include shimmy damper moments, spring restoring moments, 
inertia forces in a landing-gear structure (exclusive of the 


wheel inertia force), and spring forces arising from the 
flexibility of a landing-gear strut or of the fuselage of an air- 
plane. In general, most of these forces and moments can 
probably be considered approximately linear except shimmy 
damper moments; however, even these moments can be 
replaced as a first approximation by equivalent linear damp- 
ing moments. ‘ (See, for example, appendix A.) 

Within the scope of a linear theory, these structural forces 
and moments will depend in a linear manner on the wheel- 
center coordinates %, 6, and 7 according to expressions of 
the type 

F v ,=T 1 (D t )vz+ T,(D t )0+ T 3 (D t )y (63) 

M„= T t (p t ) v 3 + UD t )d+ T,(D t ) 7 (64) 

M„= T 7 {D t ) V3 + T s (D,)0+ T g {D,)y (65) 

where the T’s are functions of the differential operator D„ 
sometimes called transfer functions, whose specific forms 
will depend on the type of landing gear in question. 

This concludes the discussion of the forces and moments 
acting on a rolling wheel. Now these quantities will be 
utilized to set up the basic equations of motion for a rolling 
elastic wheel. 

EQUATIONS OF MOTION 
DERIVATION OF THE EQUATIONS OF MOTION 

In this section the linearized equations of motion for a 
rolling elastic wheel are set down with the aid of the equa- 
tions from the preceding sections. 

The sum of the lateral forces acting on the wheel parallel 
to the F-axis is set equal to the inertia reaction to give (see 
eqs. (20) and (27)) 

Fp,+K\(y 0 —vs—r i y)—K r 7=m a D t 3 ri3 ( 66 ) 

or, rearranging, 

F v ,+K^y 0 — (K->rr m w D t 3 )^— (K\r 3 + K y } y—0 (67) 

The first term in equation (66) is the structural force, the 
second term is the net force on the wheel resulting from tire 
elastic and inertia forces (K\=K\ tUt[lc +AK\ where A is 
given by eq. (53)), the third term is the lateral ground force 
resulting from tire tilt, and m w is the mass of the wheel (in- 
cluding the tire). For reasons previously discussed, hys- 
teresis forces and moments are not included either in this 
equation or in the following equations. 

Setting the sum of the lateral tilting moments about the 
wheel center equal to the inertia reaction gives (see eqs. 
(20), (27), (28), and (35)) 

A^,-f-F J {cxj/o+(l—Cx)’?3+[(l—cx)r 3 —c T ]7—i73} + 

[Kx(y 0 — vz~r 3 y) —K y y]r 3 —~^ Dfi=I xw D t 3 y ( 68 ) 

or 

Afa-j- (Fxr^F^yo— [Kxr 3 -t--F I c>,] D t d— 

[Kxr^+Kyn— F> s (l— cJ+Ffy+Ixw D, 3 ]y=0 (69) 
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The first term in equation (68) is the structural moment, 
the second term is tlie moment resulting from the vertical 
ground load, the third term is the moment of the ground 
forces resulting from tire lateral distortion and tilt, the 
fourth term is the gyroscopic moment resulting from the 
swiveling motion of the wheel, and I xa is the moment of 
inertia of the wheel about an axis through its center parallel 
to the X-axis. 

Setting the sum of the swiveling moments about the wheel 
center equal to zero yields the equation (see eqs. (24), (32), 
and (34)) 

M x ,+K a (p- 1 D,y o ~0)—TV DAyo-Vi—W)——- D,y=I xa D, 2 e 

(70) 

or 


wheel motion. However, before proceeding ■with these two 
matters it is convenient to digress slightly to consider the 
exact solution of the general equations for the case of steady 
yawed rolling, in order to establish several relations which 
will be useful in later sections. 

EQUATIONS FOR STEADY YAWED ROLLING 

For an untilted wheel which rolls at constant velocity at a 
constant small swivel or yaw angle, y 0 (x+h)=y 0 (x) = 
Constant, 0= Constant, and ij 3 =7=0, so that equations 
(2) (with y (1 for y t ), (16), (67), and (71) reduce, respectively, 
to the relations 

\—Vo (72) 

Vo— (L-\-K)d—l\d (73) 


d/„+ (E a v~' — rP)D l y, 1 +rv D,v 3 — (K a +I zr D, T )0— 

rr 3 )p D t y—0 (71) 

The first term in equation (70) is the structural moment, 
the second term is the net moment resulting from tire elastic 
and inertia forces exclusive of the gyroscopic moment due 
to tire lateral distortion (K a —K atlatlc +AK a where A K a is 
given by eq. (54)), the third term is the gyroscopic moment 
resulting from tire lateral distortion, and the fourth term is 
the gyroscopic moment resulting from wheel lateral tilt. 

Equations (67), (69), (71), and (16) or (19), together with 
the three auxiliary equations (63) to (65), are the basic equa- 
tions of motion for an elastic wheel. If the T-functions in 
equations (63) to (65) are known for a particular landing 
gear, these equations can be solved simultaneously to deter- 
mine the rolling behavior of the gear. 

Next the question arises as to the most profitable method 
of solution of these equations for practical landing-gear 
problems. Essentially, the choice is between exact and 
approximate solution of the equations. In the past, exact 
solutions (omitting some of the less important terms pre- 
viously mentioned) have been made only for the simplest 
case of a rigid swiveling landing gear attached to a rigid 
fuselage (refs. 2, 4, and 5). Although the exact solution 
of these equations for more complex problems does not 
appear to present any insurmountable difficulties, relatively 
complex transcendental equations may be involved, so that 
it is worthwhile to examine the possibility of finding simpler 
systematic approximations to the general equations. 

A second reason for investigating systematic approxima- 
tions to the summary theory arises in connection with the 
correlation of the summary theory with the other existing 
theories. Superficially, in its present form, the summary 
theory does not closely resemble most of the other existing 
theories. However, the approximations that are presented 
subsequently make the correlations between the different 
theories fairly easy to see. 

Subsequent sections of this report, will be concerned with 
the problem of establishing a series of systematic approxi- 
mations to the general equations and the correlation of 
these approximations with the other existing theories of 


F,,+K,y 0 =0 (74) 

M It -K a 0 = 0 (75) 

By combination of equations (72) and (73) the the lateral 
distortion is found to be 

Xo=Z,0 (76) 

By combination of equations (73) and (74) the lateral 

force on the wheel is found to be 


F^-hKtf 

The quantity l x K\, which represents the lateral force per 
unit yaw angle, is an important the characteristic called the 
cornering power or lateral guiding characteristic of the the. 
Later in this report it is found convenient to represent 
this quantity by a single symbol N, where 

N=hK x (77) 


Another property of the steady yawed rolling condition 
that is of some interest is the distance of the center of pressure 
of the lateral force behind the center of the the, which is 
sometimes called the pneumatic caster e=—M t ,/F y ,. This 
quantity, according to equations (73) to (77), is equal to 


M a K a 


(78) 


SYSTEMATIC APPROXIMATIONS TO THE 
SUMMARY THEORY 

In this section the possibilities for simplifying the preced- 
ing equations of motion are discussed, and a series of sys- 
tematic approximations to the general equations of the 
summary theory is set down. 

All but one of the equations of motion (eqs. (16) or (19), 
(63) to (65), (67), (69), and (71)) are usually simple linear 
equations and present no great difficulties. The exception is 
the kinematic equation, which was originally transcendental 
in form (eq. (16)) and was later expressed as an infinite series 
of linear terms (eq. (19a)). The most promising way to 
simplify the kinematic equation appears to be to assume that 
the series expansion in equation (19a) is a rapidly convergent 
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series in which all terms above a certain value of n can be 
neglected. The rapidity of convergence of the series and its 
significance cannot be fully determined without a knowledge 
of the particular landing-gear configuration considered. 
However, some insight into this question can be obtained by 
considering the case of purely sinusoidal oscillations of the 
form 7 /o=ri'i I ) where the quantity v y is the path frequency. 
Substitution of this expression into the infinite series in y 0 
in equation (19a) yields 


(l +Jb/n# n )2/o=(;Pi» +i]>2 a )yo 

where 

pi„ = l— W+W— • • • 

Vi<n = Uvi—Uvi J rUv y — . . . 



(79) 


(80a) 


Another form for the p’s can be obtained by substituting the 
relation y<>=e l, ' x into equation (16). The result is 


Pi„=cos vih—Lv i sin vji 
# 2 »=sin v y h-\-Lvi cos vji 


(80b) 


very rapidly. From a purely qualitative point of view the 
figure seems to indicate that, in dealing with shimmy wave 
lengths greater' than approximately 4 tire radii, two terms in 
each series are sufficient to represent fairly well the exact vari- 
ations, for wave lengths greater than approximately 6 radii 
one term in the p t series and two in the p x series are sufficient, 
and for wave lengths greater than about 20 radii one term 
in each series is sufficient. (The wave lengths cited here 
represent only order of magnitude and are not necessarily 
quantitatively significant.) To correlate these observations 
with the conditions of wave length likely to be encountered 
in practice it can be stated that the experimental data of Von 
Schlippe and Dietrich (ref. 4 or 5) and Kantrowitz (ref. 8), 
which are probably fairly typical in this respect, demonstrate 
wave lengths which are about 4 radii long at zero rolling 
velocity and which increase with increasing rolling velocity. 
Thus it appears possible that the use of only a few terms in the 
series expansion may lead to a reasonable prediction of 
shimmy characteristics for practical operating conditions. 

With the preceding considerations in mind the following 
approximations to the general wheel-motion equations were 
established. 


APPROXIMATION A 


The rate of convergence of the p series of equations (80a) can 
be tested for any given frequency by substituting numerical 
values of L, h, and v x into equations (80a) and (80b) and 
comparing the individual terms. A typical comparison is 
shown in figure 4 for the conditions £=0.8r, and A=0.5r. 
The abscissa of this plot represents the oscillation's wave 
length S—2ir/vi and the ordinate represents the p functions. 
The label p xi means that this curve represents the sum of the 
first two terms in the p Xm series, and the other labels are 
analogous. (The approximation symbols will be explained 
later.) From this figure it is seen that the series converge 
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FiauRE 4. — Variation of pi and pi with shimmy wave length. 
L=k).8r; h=Q.5r. 


As a first approximation to the general kinematic equation 
(19a) all terms for »> 3 will be neglected. This gives the 
approximate differential equation 

Dyo~\~lt D^Jo^h hoffi hd— 

r 

=v 3 +l 1 0+(r i ~^y y (81) 

This equation, together with all the general force and moment 
equations previously discussed, is referred to hereinafter as 
approximation A. 

APPROXIMATION B 

A second, less exact approximation for equation (19a) is 
obtained by setting 1^=0 for tT> 2. Thus 

Vo+h Di/o+^2 D 2 ya—Vo~{-li0 ~^— ~T 

T 

=''73+M+^r 3 — (82) 

This equation is referred to as approximation B. 

APPROXIMATION Cl 

Another, cruder approximation for the general differential 
equation (19a) is obtained by neglecting all terms in the 
series for n^> 1. This gives the differential equation 

l/o+^i Dy<i=yo-\-liQ— — - y 
r 

=’73+h®+^r 3 — - ^ y (83) 

which is referred to as approximation Cl. 
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APPROXIMATION C2 

As a slight simplification of approximation Cl, the rela- 
tively unimportant, or at least questionable, term involving 
£ may he omitted from equation (83). This gives the 
differential equation 


yo+^1 Diy Q =Tj 0 J rl l d=T] 3 -\-l l d J t-r 3 y (84) . 

which is referred to as approximation C2. 

With the aid of equations (2) and (23), equation (84) can 
be written in the more easily interpreted form 


\= — lya (85) 

or,, by using in addition equations (20), (22), and (77), as 


Ao= 


N _F v x_ NM, a 

K^K a 


( 86 ) 


Thus, in this approximation the lateral distortion of the tire 
is directly proportional to the angular distortion. 

The physical meaning of this approximation can be 
obtained by considering that equation (86) can also be 
arrived at by letting the ground-contact half-length h 
approach zero in the general differential equation (19a) (as 
was mentioned by Rotta in ref. 2), since all terms, in the 
series for »>1 and the tilt term are multiplied by h. Then 
equation (19a) (with f=0) becomes 

yo+L Dy 0 —Tj 0 -\-L6 (87) 

or, by using in addition equations (2) and (23), 


about the other tire properties. Three such assumptions 
will now be made to simplify further the equations of 
approximation C2. For the first approximation, to bo 
called approximation Dl, the simplification 

Zi=0 (91) 

is adopted. Then it follows from equation (85) that, for 
finite a, 

Xo=0 (92) 

which is the basic equation for this approximation. Thus 
for this approximation the tire is free to twist but not to 
deflect laterally. Infinite lateral stiffness is, therefore, also 
implied: 

I2x=® (93) 

For the simplest form of wheel shimmy, due to tiro 
elasticity rather than structural elasticity (to be considered 
subsequently), equation (92) does not provide accurate 
Information. For wheel shimmy due largely to structural 
elasticity rather than tire elasticity, this approximation 
may be of some value; actually most existing theories 
corresponding to this approximation have been developed 
for the primary purpose of considering the influence of 
structural elasticity on wheel shimmy. 

APPROXIMATION D2 

- As a second simplification of approximation C2, tho 
assumption 

Zj=co (94) 


X=-Aa 

Also, equation (77) for the yawed rolling becomes 


can be adopted. The corresponding theory is designated 
as approximation D2. From equation (85) it is evident that 
this approximation implies that, for finite 


N—K\L 


(89) 


a=0 


(96) 


and the combination of equations (88) and (89) gives 


which in turn implies that 


(90) 

Equation (90) is the same as equation (86) for any given 
combination of N and K\. Thus, when written in the form 
of equation (86), approximation C2 formally corresponds to 
the assumption of !i=0. 

Reliable qualitative results should be expected from 
approximation C2 only when the neglected quantity h is small 
with respect to the characteristic length of the rolling motion 
in question (for example, the wave length S of a sinusoidal 
oscillation). Fortunately, this condition is at least some- 
times satisfied for practical rolling conditions. 

APPROXIMATION Dl 

Before considering the next approximation it should be 
remembered that all of the terms neglected in the preceding 
approximations were multiplied by the tire ground-contact 
half-length h; thus these approximations implied the assump- 
tion of progressively smaller ground-contact length or 
progressively larger wave length. In order to simplify the 
equations further, it is necessary to make some assumptions 


N= o° 

K a =eN=<*> (96) 

Thus for approximation D2 the tire is considered to bo 
torsionally rigid but laterally flexible. 

APPROXIMATION D3 

A third simplification of approximation C2 can be obtained 
by keeping the quantity li finite but considering the tiro to 
have both infinite lateral stiffness and infinite torsional 
stiffness, or 

E x =K a =N= co (97) 

This approximation, which is designated as approximation 
D3, thus represents the case of a rigid tire and consequently 
also implies that a=X o =0. (Formally, approximation D3 
can also be interpreted as the limiting subcase of approxima- 
tion Dl where N= co or as the limiting subcase of approxima- 
tion D2 where Kx= ® . However, it should not be con- 
cluded that approximation D3 is necessarily inferior to these 
other two approximations.) 

A choice of seven simplified approximations based on the 
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summary theory is now available. The problem remains of 
determining which, if any, of these approximations is the 
simplest one which can be used for any particular tire- 
motion problem. While it is not yet possible to give a 
completely satisfactory answer to this question, some insight 
into the answer can be gained, by comparing the various 
approximations with the previously published tire-motion 
theories, which are at least partly successful and most of 
which are closely related to these approximations. 

CLASSIFICATION AND EVALUATION OF EXISTING 

THEORIES 

In this section the major previously published theories of 
wheel motion are briefly reviewed, evaluated, and, wherever 
possible, correlated with the preceding summary theory of 
this report and its approximations. Each of the major 
existing theories is first considered individually and then an 
abbreviated overall summary classification is presented in 
tabular form. 

INDIVIDUAL REVIEW AND EVALUATION OP EXISTING 
THEORIES 

VON SCHUPPE AND DIETRICH 

The tire-motion theory of Von Schlippe and Dietrich 
(refs. 3 to 5), of course, corresponds directly to the-summary- 
theory of this report, since the summary theory was taken 
from their theory with only minor modifications. These 
modifications included a more detailed consideration of 
some of the influences of lateral tilt and of tire inertia forces 
and moments. It should be noted, however, that the Von 
Schlippe-Dietrich theory is more advanced than the sum- 
mary theory of this report in that it partly takes into account 
the width of the ground-contact area. However, as was 
previously noted, this factor is probably not of great practical 
importance. 

ROTTA 

Rotta’s tire-motion theory (ref. 2) corresponds to the 
summary theory of this report because it also is based on the 
Von Schlippe-Dietrich theory. Rotta’s theory represents a 
slight extension of the last theory to take into account more 
adequately most of the effects of tire tilt and the width of the 
ground-contact area. No inertia forces due to tire lateral 
distortion or centrifugal forces are discussed. 

BOORCIER DE CARBON ADVANCED THEORY 

Bourcier do Carbon (ref. 6) has developed two closely 
related theories of tire motion which are similar to approxi- 
mations B and C2. The first of these will be referred to as 
the Bourcier de Carbon advanced theory and the second as 
the Bourcier de Carbon elementary theory. 

Bourcier de Carbon’s advanced theory uses five basic tire 
properties which are correlated with those of the present 
report through the following relations: 

D=~ (98a) 



(98c) 


(98d) 

w 

•p li 2 ( L-\-K ) 

hK a KJi(2L+h) 

(98e) 


Equations (98) were obtained by comparing this theory 
with the corresponding approximation B. The symbols of 
Bourcier de Carbon are overscored and do not necessarily 
bear any relation to other symbols in this report designated 
by the same letters. Although the symbols D, T, S, and 
e bear a simple relation to the symbols used in the deriva- 
tions of the present report, the symbol R bears a more com- 
plicated relation which is worth some detailed consideration. 

Bourcier de Carbon defined the tire property R as follows: 
If an untilted wheel is rolled forward while exposed to a 
constant turning moment about a vertical axis and with no 
side force, it will move in a circular path of a definite radius; 
R is defined as the reciprocal of the product of the turning 
moment and the path radius. Unfortunately, however, this 
constant-moment circle-rolling experiment is not easily 
performed. Therefore, equation (98e), which expresses R 
in terms of the more easily measured fundamental quantities 
L, h, and K a , is of importance for the use of the Bourcier de 
Carbon advanced theory. 

In treating the subject of tilt, Bourcier de Carbon omits 
many of the details considered in this report. For example, 
he implicitly assumes that K 7 =c\—c 7 =£—0 and that the 
inclination angle k is small (taking cos k~ 1). ' However, 
these omitted tilt terms may be as important as the terms 
considered (as will be shown later); therefore, Bourcier de 
Carbon’s considerations of tilt are incomplete. 

It should be noted that in reference 6 certain misconcep- 
tions occur in the parts of the paper that deal with com- 
parisons between theory and experiment. In particular, 
some of the experimental data quoted by Bourcier de 
Carbon from reference 3 of the present report appears to 
be either misquoted or misinterpreted. Consequently, 
Bourcier de Carbon’s conclusion that the experimental data 
of reference 3 provide a remarkable check of his theory is 
not completely justified; actually these experimental data 
provide only a fair indirect check of the theory. 

GREIDANTJS 

Another theory similar to approximation B, except for 
the influence of tilt, is that of Greidanus (ref. 7). Greidanus 
considers the influence of tilt in much greater detail than 
does Bourcier de Carbon; however, he also fails to consider 
the force term proportional to K y , and thus his results also 
do not fully describe the influence of tilt. 

In addition, Greidanus’ kinematic equation differs from 
equation (82) for approximation B in that he has introduced 
a slightly different term associated with tilting of the tire. 
In the present terminology Greidanus’ equation reads 



y 0 +li Dy 0 +1 a D 2 yo=’?o+li^ — 


(98b) 


(99) 
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The difference between the two equations lies in the co- 
efficient of 7. For approximation B (eq. (82)) the coefficient 
is 


£Lh 

r 


( 100 ) 


and for Greidanus’ equation (after substituting for h from 
eq. (19a)), 

( i+ 2 ) h (101) 


If £ is set equal to (^L+ jlL, 


the two coefficients are iden- 


tical; thus Greidanus’ kinematic equation can be considered 
to be a particular case of the corresponding equation of 
approximation B. 

No subsequent detailed discussion of Greidanus’ theory is 
included in this report because no complete translation of 
reference 7 is available. 


these versions is governed by the equation 

a+TD,a=-^f=~^ (102) 

or 

Tkv IPyM Dy 0 +y 0 =vo+ Tl x v Dd+he (103) 

where T is a time-lag constant. . This theory corresponds to 
a generalization of approximation C2 (with pneumatic caster 
neglected, that is, 6=0), since for T= 0 equation (102) is 
identical with the basic equation for approximation C2. 
However, for 7V0 this theory is not directly compatible 
with the s umm ary theory and its approximations. 

Moreland uses the following reasoning to establish this 
equation: First, it is known that in steady yawed rolling a 
yaw angle a. is developed as a consequence of the application 
of a lateral force F„\ according to the relation 

(104) 


BOURCIER DE CARBON ELEMENTARY THEORY 

Bourcier de Carbon's elementary theory corresponds to 
approximation C2 of this report except for the minor short- 
comings which were discussed in connection with the Bourcier 
de Carbon advanced theory. The only difference in Bourcier 
de Carbon’s two theories is that the coefficient R is taken as 
infinity in the elementary theory but is finite (see ecp_ (98e)) 
for the advanced theory. The infinite value for R corre- 
sponds to the assumption Z 2 =0, which was previously made 
in passing from approximation B to approximation C2 
Compare eqs. (82) and (84)). The physical significance of 
R—m is obvious from equation (98). It means that h= 0. 

MELZER 

The Melzer theory of tire motion (ref. 10) is also similar 
to approximation C2 except for details of the tilting process. 
Otherwise Melzer’s kinematic equation is identical with the 
kinematic equation of approximation C2 and of Bourcier de 
Carbon’s elementary theory. However, Melzer’s theory 
differs in that it treats the moment due to tire twist as pro- 
portional to the swivel angle —6 rather than to the tire 
twist angle Dy 0 —6. This assumption would appear justified 
only if Zty o <C0, which is not true in general. It is interesting 
to note that for the simplest case of wheel shimmy (see sec- 
tion entitled “Application to Wheel-Shimmy Problems — 
Case I”) the Melzer approximation leads to one of the same 
stability boundaries and to the same limiting high-speed 
shimmy frequency as the more nearly correct approximation 
that included the term in Dy 0 . This restricted agreement, 
however, hardly justifies the use of Melzer’s approximation, 
since predictions made by means of the two approximations 
differ with respect to divergence of the shimmy oscillations 
and with respect to another stability boundary. Moreover, 
for simple problems the Melzer approximation is not sig- 
nificantly easier to solve than the more nearly correct form 
including the Dy 0 term. 

MORELAND ADVANCED THEORY 

Moreland has proposed three versions of a tire-motion 
theory in references 11 and 12. The most advanced of 


which is the basic equation for approximation C2. How- 
ever, for the dynamic rolling case this equilibrium yaw angle 
obviously cannot be established immediately upon applica- 
tion of a given side force; rather, a finite amount of time will 
be required for the equilibrium yaw angle to develop. 
Moreland has attempted to take this finite time lag into 
account by modifying equation (104) to the new form of 
equation (102). In the latter equation the constant T is a 
measure of the- time lag of the yaw angle behind the applied 
force Fj\. 

This time-lag term introduced by Moreland does not 
correspond exactly to any of the terms in the summary 
theory, and to this extent Moreland’s advanced theory is 
apparently incompatible with the summary theory. How- 
ever, a partial reconciliation of the two theories can bo 
obtained through the following considerations of hysteresis 
effects as applied to" approximation C2 : According to equa- 
tion (58) the lateral ground force, if tilt and inertia forces 
are neglected, is given by the equation 

X 0 +2 7 xAXo=^ (105) 

and the kinematic equation for approximation C2 is (eq, (85)) 

Ao= — 

Combining equations (105) and (85) to eliminate Ao and 
substituting Kyk=N yields the equation 

a+T x r> ( a=-^ (106) 

Equation (106) is formally identical with Moreland’s basic 
equation (102) if the hysteresis time constant T\ is con- 
sidered equivalent to Moreland’s time constant T. The 
important points to be noted here are: (1)' according to both 
views, the tire twist a lags behind the applied lateral force 
F t \, and (2) Moreland adopts the lateral-force equation 
F v \—KxSq which implies that the lateral force and lateral 
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deformation are in phase (although lateral force and twist 
are not in phase), whereas according to approximation C2 
(see eq. (105)) the lateral force lags behind the lateral de- 
formation as a consequence of the kinematic relation X= — l^a. 
In regard to the second point, since there are apparently 
no pertinent experimental data available, it is not possible 
to conclude w r hich point of view, if either, is correct. 

Why Moreland’s time-constant concept has not been in- 
corporated directly into the derivations of the summary 
theory should be clear from the preceding discussion: It is 
not certain whether Moreland’s time-constant terms are 
really independent of the terms already contained in the 
summary theory or whether they are, rather, another way 
of looking at some terms which are already included in the 
summary theory. More specifically, Moreland’s analysis 
does not include inertia forces and moments due to tire 
lateral distortion, hysteresis forces and moments, or the 
higher l n terms (lj, l 3 , . . .), and for certain conditions any 
of these factors could be interpreted as a time-constant 
effect. In view of these factors and in view of the lack of 
pertinent experimental data, a completely satisfactory eval- 
uation of the relative values of the summary theory and 
Moreland’s advanced theory cannot be made in the present 
report. 

MORELAND INTERMEDIATE THEORY 

As a simpler approximation for his advanced theory, 
Moreland has implied (ref. 11) that the influence of the time- 
lag term in the basic equation for his advanced theory (eq. 
(102)) can be approximated for the usual range of shimmy 
frequencies by using the simpler kinematic equation 

40Z I a=— M (107) 

Inasmuch as approximation C2 has the kinematic equation 
(85) 

li<x= — Xo 

and approximation D2 has the kinematic equation (95), 
which could be written in the form 

— Xo 

it follows from a comparison of these last three equations' 
that Moreland’s intermediate theory falls somewhere be- 
tween approximations C2 and D2. Since Moreland has not 
offered any concrete justification for this approximation, 
further detailed discussion does not appear warranted here. 

MORELAND ELEMENTARY THEORY 

Moreland’s most elementary theory corresponds directly 
to approximation D3, the case of a completely rigid tire, ex- 
cept that it, like Moreland’s other two theories, does not 
take into account the pneumatic caster (e=0). 

TEMPLE ELEMENTARY THEORY 

Temple has proposed an elementary theory for the motion 
of tires which is identical with approximation Dl (ref. 13). 
Temple has chosen the most general form of this approxima- 
tion in that he has considered both the tire torsional stiffness 
K a (indirectly interpreted as an increase in trail) and the 
cornering power N. 


1.59 

This theory was developed before experimental evidence 
pointing to the need for more detailed considerations of tire 
lateral stiffness was available. Subsequently, Temple has 
indicated a need for more refined considerations of the tire 
(ref. 25) and has developed independently a theory (unpub- 
lished, but partly described in ref. 21) similar to the theory 
of Yon Schlippe and Dietrich. 

MAXBR 

Maier (ref. 14) has proposed a simplified theory similar to 
approximation Dl, with the difference that he makes the 
added assumption that the tire torsional stiffness K a is zero. 
This theory, too, was developed before there existed much 
experimental evidence pointing to the need for more refined 
considerations for shimmy behavior. 

TAYLOR 

Taylor (ref. 15), in -a brief paper, suggested another tire 
motion theory which corresponds to approximation D2 
except that details of the tilt process are omitted. 

KANTROW1TZ AND WYLIE 

The preceding theories for tire motion, which include 
most of the known theories, may all be considered as closely 
related to the summary theory of this report. However, 
two other well-known theories, one by Kantrowitz (ref. 8) 
and one by Wylie (ref. 9), apparently cannot be derived 
from the summary theory and thus cannot be accurately 
classified here with respect to the other theories. They 
possess some of the merits of approximation B but in other 
respects are less adequate than approximation 02. To 
point out the deficiencies of these two theories it is sufficient 
to consider two simple cases of tire motion as follows: 

The first case to be considered is the steady straight-line 
motion of a nonswiveling, untilted, rolling wheel which is 
not yawed with respect to its direction of motion (a=0), 
which is inclined by an angle 9 (equal to the swivel angle) 
to the reference Y-axis, and which has no lateral forces or 
moments acting on the wheel. (See sketch 4.) Obviously, 



Y 


Sketch 4 

for this case there will be no lateral distortion of the tire, or 

Xo=0 

On the other hand, Kantrowitz’ basic equation, which is 
Xo"i“D DXo— "LB — 12 T)9 (108) 


400104 — 0 ! 


•12 



REPORT 1299 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


160 

gives for this steady Tiny awed case (with D\=D9= 0) 

\ —Le 

This equation is obviously incorrect, since it implies that the 
lateral distortion of a straight rolling wheel, which actually 
must be zero, depends on the choice of the coordinate axes 
to which 0 is referred. Only for the special case of a wheel 
running along the reference axis (that is, for 6=0 ) is 
Kantrowitz’ theory correct in this respect, and in an actual 
shimmy problem this condition is possible only for the case 
of zero trail; hence Kantrowitz’ theory cannot necessarily 
be expected to give reliable results for trails different from 
zero. Thus Kantrowitz’ theory is of at least doubtful value 
for practical shimmy calculations. 

In order to evaluate the Wylie theory, consider the case of 
steady untilted yawed rolling of a wheel moving parallel to 
the X-axis. (See sketch 5.) Obviously, the lateral distor- 
tion of the tire \ will depend only on the swivel .angle 9 



Y 


Sketch 5 

(9= a) and not at all on the absolute lateral displacement of 
the wheel y 0 - On the other hand the basic equation of 
Wylie, which in the present terminology is 

y 0 +L Dy 0 =L9—ls D9 (109) 

gives for this steady case (where Dy a =D9=0) the relation 
y^=Jj9 or, by using equation (2), 

—L9 — ij 0 

This equation states the obviously incorrect conclusion that 
the tire distortion is dependent on ij 0 or, in other words, that 
it depends on the choice of the coordinate axes. Thus, only 
for the special case 7?o=0 is Wylie’s theory plausible in this 
respect, and j / 0 =0 implies that the reference axis must pass 
through the path of the wheel. Since this condition is 
satisfied in to actual shimmy motion only for the special 
case of zero trail, Wylie’s theory, like Kantrowitz’, can be 
fully valid only for zero trail and, consequently, this theory 
is also of doubtful value for practical shimmy calculations. 

It might he noted that the preceding difficulty concerning 
Wylie’s theory could be removed by the logical procedure 
of adding the term ij 0 to the right-hand side of Wylie’s 
equation (109) to give the new basic equation 


Another questionable point in these two theories is that 
the Kantrowitz theory, as -previously noted, predicts that 
the lateral distortion in yawed rolling parallel to the X- 
axis (see sketch 5) is 

Xo =L9 

and so does the Wylie theory if the reference axis is chosen 
to give 7? 0 =0. On the other hand, the summary theory has 
for this yawed case the relation 

Xo= (L-\-K)9 

(See eq. (85).) The difference arises from the fact that 
Kantrowitz and Wylie did not explicitly consider the ground- 
contact length 2 h in their derivation. 

OTHER THEORIES 

In addition to the just discussed theoretical papers dealing 
particularly with the subject of landing-gear shimmy, a 
number of relevant papers exist which are either largely of 
historical interest, which deal, particularly with automobile 
shimmy problems, which deal only briefly with landing- 
gear shimmy problems, which deal with other tire-motion 
problems such as yawed rolling and veering-off or ground 
looping, or which deal with the determination of tire stiff- 
ness parameters. Although these papers are of some interest, 
they do not appear to contain any important contributions 
which are not contained in the theories just reviewed. The 
reader is referred to reference 1 for a substantially complete 
listing and brief discussion of most of the papers in this 
class. 

Of particular historical interest among the investigations 
not considered here in detail are the work of Broulhiet (ref. 
23) and the work of Fromm (discussed in ref. 22). These 
two investigators independently were apparently the first 
to recognize the importance of lateral distortion and corner- 
ing power of tires in the wheel shimmy problem. Taking 
these factors into account, both authors developed tire- 
motion theories whose kinematic relations correspond to 
that of approximation C2 of the present report. 

TABULAR CLASSIFICATION OF EXISTING THEORIES 

In order to permit easier visualization and comparison 
of the merits of the theories discussed, the major assump- 
tions of the various theories of tire motion are collected in 
table I. This table lists the nature of the assumptions 
made in regard to the primary tire parameters N, K\, K a , 
6, and l K for each of the theories discussed. 

APPLICATION TO WHEEL-SHIMMY PROBLEMS 

In the preceding sections of this report a set of basic 
differential equations for the motion of an elastic wheel has 
been derived and compared with the corresponding equations 
of. most of the previously existing theories. Those com- 
parisons have indicated that, from a mostly qualitative point 
of view, the summary theory of this report and the system- 
atic approximations to it incorporate the major merits of the 


ya-j-L Dy a =L6—l? DO-\-Vo 


( 110 ) 
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TABLE I.— PRIMARY ASSUMPTIONS FOR THE VARIOUS THEORIES OF TIRE MOTION i 


Theory 

N 

Kx 

Km 

MM 

a 

a 

a 

I«(n>3) 

Remarks 

Summary theory 

F 

F 

F 

F 

F 

F 

F 

wm 


Approximation A 

F 

F 

F 

F 

F 

F 

F 



Approximation B- 

F 

F 

F 

F 

F 

F 

0 

0 


Approximations 01 and C2 

F 

F 

F 

F 

F 

0 


0 


Approximation Dl— 

F 

09 

F 

F 

0 

KM 


0 

Assumes laterally rigid tire 

Approximation D2 

CO 

F 

00 

F 

CO 

IBS 


0 

Assumes torsionally rigid tire 

Approximation D3 

CO 

09 

CO 

F 

F 

■9 

0 

0 

Assumes laterally and torsionally rigid tire 

Von Sohllppo-DIotrlch and Eotta 

F 

' F 

F 

F 

F 

n 

n 

F 


Bourdor do Carbon advanced 

F 

F 

F 

F 

F 


IBS 

0 

- 

Oroldanus.. 

F 

F 

F 

F 

F 


•19 

0 


Boarder de Carbon elementary. 

F 

F 

F 

F 

F 


0 

0 


Molicr 

F 

F 

F 

F 

F 


0 / 

0 

Assumes tire twist angle=swlvel angle 

Moroland advanced _ 

F 

F 

0 

0 

F 


0 


Introduces tlmo-eomtant term 

Moreland Intermediate _ 

F 

F 

0 


F 




Implies extremely largo Ii value 

Moroland elementary ... 

CO 

CD 

CO 


F 


0 


Assumes laterally and toratonally rigid tiro 

Temple elementary.- 

F 

00 

F 




0 


Assumes laterally rigid tire 

Major 

F 

CO 

0 


0 

0 

0 


Assumes laterally rigid tire 

Taylor 

00 

F 

CO 


CO 

0 


0 

Assumes torsionally rigid tire 


F 

F 

F 


F 

0 


0 

(“For trail not equal to xero both of these 

Wyllo 

F 

F 

F 


F 

0 

■9 

0 

I theories can lead to erroneous conclusions 


i The symbol F Indicates o Unite number. 


existing theories of tire motion and avoid some of their dis- 
advantages. However, it still remains to investigate how 
best to apply the theory to specific landing-gear problems, 
to investigate the question of the absolute or quantitative 
accuracy of the summary theory and of the other theories) 
and, if the summary theory be found satisfactory, to es- 
tablish the simplest systematic approximation to it which 
will give reliable information regarding any particular prob- 
lem in tire motion. The best way to accomplish these 
various aims appears to be through a discussion of the 
shimmy of several particular landing-gear configurations. 
In this section three particular landing-gear configurations 
are discussed which range in complexity from the simplest 
case of a rigid swiveling landing gear to the most general 
case of a gear of arbitrary complexity. 

DESCRIPTION OF PARTICULAR CASES CONSIDERED 

The first landing-gear configuration considered, which is 
designated case I, is illustrated in figure 5. It consists of a 
rigid landing gear whose only degree of freedom other than 
tiro distortion is rotation of the wheel about an inclined 
swivel axis, which may bo opposed by a linear spring or 
damper. This particular configuration is chosen because 
most of the existing experimental data have been obtained 
for such a configuration. Thus, this configuration makes it 
possible to discuss and evaluate the summary theory, its 
systematic approximations, and the existing" theories with 
respect to agreement with experiment in regard to the various 
important characteristics of a shimmy motion, such as 
stability boundaries, shimmy frequency, and divergence. 

The second landing-gear configuration, case II, is an 
untilted landing gear possessing two degrees of freedom aside 
from tire distortion. This landing-gear configuration, which 
is illustrated in figure 6, consists of a wheel free to swivel but 
not to tilt, which turns about a rigid vertical swivel axis 
attached by a spring k to the supporting structure. (This 
spring is an idealized representation of the lateral flexibility 
of an actual landing-gear strut.) This configuration is dis- 
cussed for two purposes: (1) to illustrate the effect of struc- 
tural elasticity on wheel shimmy behavior and (2) to provide 



Figube 5. — Configuration of landing gear for case I. 

an example which is better suited than case I for bringing out 
the relative merits of several of the systematic approxima- 
tion theories for a case involving structural flexibility. 

The third landing-gear configuration considered is a 
modification of the gear of case H. In case H the landing 
gear was considered to be connected to its supporting struc- 
ture by a single spring; in case III this single spring is 
replaced by a more complex structure described by some 
transfer function. This case is chosen mainly to demon- 
strate the applicatiqn of the theory to complex problems for 
which the transfer-function concept may be of value. 
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Figure G. — Configuration of landing gear for case II. 


CASE I 

GENERAL DERIVATION 

In this section the basic equation of motion will be derived 
according to the summary theory for the special case of an 
inclined, rigid, swiveling landing gear (case I), which is 
illustrated in figure 5. This equation of motion could be 
obtained by making use of the previously derived equations 
of motion for the completely general case; however, it can 
more easily be derived here in a slightly different form for 
this particular problem. 

The geometric quantities which enter the discussion of this 
particular landing gear are indicated in figure 5. This gear 
has a swivel axis lying in the XT-plane and is inclined 
forward from the vertical Z-axis by a constant angle k. The 
perpeudicular distance a between the ground-contact center 
point 0 and the swivel axis is called the trail. The swivel 
axis is assumed to move with constant velocity v along the 
A'-axis without lateral motion from the XT-plane. 

[Rotation of the wheel structure about the inclined swivel 
axis by an amount p results in a component 6 of angular 
rotation about the vertical axis of magnitude 

8 = p COS K (111) 

a component y of rotation about the X-axis (tilt) of magni- 
tude 

y= — p sin k (112) 

and a lateral deflection % of magnitude 

rto=—ap (113) 

where all angles except k are considered small. 

The sum of all moments about the swivel axis must equal 
the inerfia reaction J#. D l 2 p=D??D 2 p, where 1$ is the moment 
of inertia of the wheel structure (including the wheel) about 
the swivel axis. ' The moments about the swivel axis are 
assumed to consist of the moments resulting from the 
previously discussed forces and moments that arise from 
the distortion and ground loads plus the moments applied 
to the wheel by the supporting structure, which are assumed 


to consist of a restoring spring of moment pp and a linear 
damper of moment g D,p=gv Dp, where p and g are con- 
stants. Thus, summation of the moments about the swivel 
axis gives the differential equation 

— [K\ (y Q —vo) —K y y]a—F,[<xy 0 + (1 — Ca)ij 0 — <yy] sin k+ 

K a (Dy 0 —d) cos K—r^D(y 0 —ijo) cos k— pp— gvDp=I^,v > D i p 

(114) 

where the first term is the total ground force due .to tire 
lateral distortion and tilt (see eqs. (20) and (27)) times its 
moment arm a; the second term is the vertical force times 
its moment-producing fraction sin k times its moment arm 
(see eq. (28)); the third term is the moment about the 
Z-axis due to the twist (see eq. (24)) corrected by cos \p for 
the component about the swivel axis; the remaining terms 
on the left-hand side represent the gyroscopic torque due 
to lateral the distortion (see eq. (32)), the spring restoring 
moment, and the linear damper moment. Now by making 
use of equations (111) to (113), equation (114) can bo 
written in the form 

Ai D 2 p-\-A 2 Dp-\-A z p-\-B\ Dy o -\-B 2 y o =0 (116a) 

where 

A 1 =Itt? 

A 2 — an? cos K+gv 

(116b) 

A 3 ~a 2 K\-\-K a cos 2 K-f-p-hp* *■ 

v 

Bx——K a cos K+TP 2 COS K 
B 2 =aK\+c\F z sin k 

and 

p z —aKj sin K—aF, sin k-\-oc\F, sin k-\-c v F, sin 2 k (116c) 

The general relation between p and y 0 for this case is found 
by substituting for %, y, and 6, according to equations (111) 
to (113) in the general kinematic equation (19a). Thus 

2/o+S lnD n y 0 =—ap+l 1 p cos k+^^P sin k 
n-1 T 

or, abbreviating, 

tr= 1+^ tan k (116) 

and rearranging, 

(trl-! cos k— a)P=y 0 +'^2 InD^jo^^h, In D*Vo (117) 

noo 

since lo=l. Differentiating this result gives 
(crli cos k—o)DP=F>~', l„ D n+i y 0 

7J = 0 

which can also be written as 

(cli cos K—a)DP=y~', D - < D K u n (118) 

n°l 

and, similarly, 

(cZi cos k — u)DV = 2H Z h _ 2 D'yo (119) 
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Substituting these relations into equation (115) and multiply- 
ing through by <sl x cos K—a- gives 

A x ln -2 D n V o+-42 S ^n— i D n yo+A i 2D Z„ D*y o+ 

n=2 n-l ' n=0 

B x (d x cos K—a)Dy Q +B a (d x cos K—a)y 0 = 0 

Finally, after adding all terms of like order, substituting 
N=l x K\ (eq. (77)), substituting for some of the A’s, and 
using equation (116), there results the equation 


S Fn D*y 0 =0 ( 120 a) 

n=0 

where 


F 0 =<raN cos K+K a cos 2 k+p+p k +u^ 


F x =a 2 N+a'K a cos K+pl x +pJ, x +gv+d x rF cos 2 k 
F a =A i +AjL x +A a l a 


( 120 b) 


F n = A x l n ^ x + AJ„- X + A 3 l n (n+ 2 ) J 


and 


u t —c\F,(d x cos K—a ) sin k 
a'=a+(l— <r)h cos ic=a— ^^sin k 


( 120 c) 


SYSTEMATIC APPROXIMATIONS 

Equations ( 120 ) provide the general differential equation 
of free motion for case I. The corresponding equations for 
the systematic approximations A to D3 are obtained as 
follows. 

Approximation A. — The basic equation for approximation 
A is obtained by setting Z„ equal to 0 for 3 in equations 
(120). The resulting equation can be stated in the following 
form: 


E> D%+E x D'yo+Et IByo+Es TPyo+E^ Dy 0 +E 6 y „=0 ( 121 a) 
where 


Eq — I^Flz 

E x =I^v 2 l J +(aTV i cos K+gv)l z 

E % — Wk + (art? cos K+gv)h+ (a 2 K+ 

K a COS 2 K + p + p t )k 

E 3 =I}V t + (art? cos k+<7»)Z 1 +(u 2 Zx+ 

K a COS 2 K + p + pt)h 

Ei=a 2 N+a'K a cos k+ pl x + pj x +gv+ dirt? cos 2 k 

Eq = aa N cos K+K a cos 2 k+ p+ p c +u x 

and 


( 121 b) 


Approximation B. — In order to obtain approximation B, 
set 1-3 equal to 0 in equations (121b). This gives the differ- 
ential equation 

Eq TFijq-j-Ei IByo+Ei D i yQ+JE } DyQ+E x y 0 =0 ( 122 a) 

where 

Eq—I+FI^ 

E l =I^Fli+(arF cos K+gv)^ 


Ei= I+F+ (arF cos K+gv)h + (a 2 K\+K a cos 2 /c+p+ pl)h 


H 122 h) 


E 3 =a 2 N-\-a'K a cos k+ pl x + p++gv+ d x rF cos 2 k 
E±= aaN cos K+K a cos 2 K+p+p t +u c 


(The corresponding equation with inclusion of hysteresis 
effects is listed in appendix B.) 

Approximation Cl. — The equation for approximation Cl is 
obtained by setting U equal to 0 in the equation for approxi- 
mation B. The resulting differential equation is 

Eq IFya+Ei Dhjo+Et Dijq+E^Jq^O (123a) 


where 


Eq=I+FIi 


Ei=Ii,F+(arF cos K+gv)l x 

Et—a^N+a'Ka cos k+ pl x + pji x +gv+ d x rF cos 2 k 


(123b) 


E^-craN cos K+K a cos 2 K+p+p x +u K . 


Approximation C 2 . — The equation for approximation C2 is 
obtained by setting £ equal to 0 in the equation for approxi- 
mation Cl. The resulting differential equation is 

K IFyo+Ei D^q+Ei Zh/o 4" £ 32 / 0=6 (124a) 

where 


E 0 =/*a 2 Zi 

Ei = I 4 F + (an? cos K+gv)l x 

E 2 =a 2 N+aK a cos k+ pl x + pJ, x +gv+l X Tir cos 2 k 

Ei=aN COS K+K a COS 2 K+p+pt+U t i 

and 

u kX =c\F,(I x cos K—a) sin k 


(124b) 


(124c) 


(The corresponding equation, with hysteresis effects in- 
cluded, is listed in appendix B.) 

Approximation Dl. — The differential equation for approxi- 
mation Dl is obtained by setting l x equal to 0 in the equation 
for approximation C 2 . The result is the differential equation 


Eq I+ya+Ei Dyo+Eiy 0 =0 (125a) 

where 

Eq=I*F 

E 1 =a 2 N+aK a cos K+gv 
EQ=aN cos K+Ka cos 2 K+p+w x J 


p,= (aK y —aF t +acy.F t +c 7 F, sin k) sin k 
u t =c\F t (d x cos K—a) sin k 


( 121 c) 


^ (125b) 
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and 

/ w t ={aKy—aF I -\-c y F t sin k) sin k (125c) 

Approximation D2. — The differential equation for approxi- 
mation D2 is obtained by dividing all terms in equations 
(124) by li and then setting l t equal to co and using the 
relations Nfli=K\ and K a !h=eK\. The resulting equation is 

Ed IFy a -\-Ei IPy 0 +E Dy 0 -\- Eyo ■= 0 (126a) 

where 

Eo=Itx? 

Ei=arv* cos K-\-gv 

Ei = a?K\ + cuK\ cos K+p+Pr+ro 2 cos 1 k 

E 3 =aK\ cos K+eK\ cos 2 /c-f c\F, sin k cos k 

Approximation D3. — The differential equation for approxi- 
mation D3 is obtained by dividing the E’s for approximation 
Dl by N and then setting N equal to <*> . The resulting 
differential equation is 

a Lh/o+y 0 cos k=0 (127) 

PREVIOUSLY PUBLISHED THEORIES 

In the preceding section the differential equations are 
set down for case I according to the summary theory and 
the systematic approximations thereto, complete with all 
pertinent details, including a number of second-order terms 
so as to enable the direct application of these equations to 
actual problems. In the present section the differential 
equations for the previously published theories are listed 
for the purpose of making clear the differences in the basic 
structure of the differential equations resulting from the 
application of these theories to case I. In order to avoid 
obscuring the more important differences between the 
equations of the various theories, all terms are omitted 
whose inclusion in any shimmy theory should be completely 
straight-forward (such as the spring restoring-moment 
coefficient p and the damping coefficient g ) and also all 
inclination effects. Although the latter effects are not 
necessarily negligible, they do appear to be of second-order 
importance and their omission here should not alter the 
basic structure of these equations. With these omissions — 
that is, for p=p=/c=0— the differential equations for case I 
according to the previously published theories are as follows 
in the terminology of the present report: 

Von Schlippe-Dietrich and Rotta theories. — The basic equa- 
tion of motion for this case according to the Von Schlippe- 
Dietrich theory, after neglect of the effects of tire width, 
corresponding to equation (120a) for the summary theory, is 

T;F a 'D”y { 0 =0 ■ (128a) 

n-0 



where 

F 0 '=aN+K a 
Fi =a 3 N+aK a +g 1 v 
E n ' =IiFl n -iFgivl n -i-\-(a 2 K\ J rK a )l n -\- 

[ja! dh ' K ^~ a )\ 

(n=2, 4, ... ) 

F {a 3 K\+K^)l » — 

(n= 3, 5, ... ^ 

The symbol pi is a hysteresis damping coefficient defined 
by the equation 

(128c) 

where di and c/ 2 are hysteresis constants. Aside from tho 
omission of some inclination effects and other terms, as 
discussed previously, the only differences between equations 
(128) and equations (120) for the s umm ary theory lies in 
the inclusion of the hysteresis term involving g x and in the 
addition of the terms in brackets in equations (128b). 
Rotta’s corresponding equations, after neglect of tire width 
effects, would bo the same as equations (128) except that 
Rotta omitted the hysteresis term. As was previously 
noted, the Von Schlippe-Dietrich theory and Rotta’s theory 
differ only slightly in their respective considerations of tho 
influences of tilt and tire width, neither of which effects are 
considered here. 

Bourcier de Carbon advanced theory. — The Bourcier do 
Carbon advanced theory leads to the fourth-order differential 
equation 

Eq D l y(,~\~Ei IFyo^Ei Dhy 0 -\-E Dj/o'h-Ed/o— 0 (120a) 

where 

E 0 =I+v‘li 

Ei=IfFLi 

E 3 =I^+(a i K x +K a )h 
Ei — (zW -f- (il\ a 
Ei=aN+K a 

Bourcier de Carbon elementary theory. — Tho coefficients 
for the Bourcier de Carbon elementary theory are obtained 
by setting U equal to 0 in equations (129) for the advanced 
theory. The resulting third-order differential equation is 

Ed TFyD~\~Ei D^q-^Ei Dyo-\-Eiy 0 =0 (130a) 
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wliore 


E=I^k 

E=E# 

E i =a 2 N+aK a 

Ez=aN+K a 


(130b) 


Melzer. — Melzer’s theory gives the third-order differential 
equation 


where 


2?o IE/oEEi E^y^E Dyo+Eyo — 0 

l?o=/ t ttr 1 Z 1 

Ei=7#o» 

E i =d‘N+liK a f 


(131a) 


(131b) 


where 

E=Irf •) 

E = a W+ aK a >- (135b) 

E=aN+K a J 


Maier. — The second-order differential equation of Maier’s 
theory is, 

Ei EE)(i~\~ Ei Dy<f\~ EUo—Q (136a) 

where 

E=o?N [- (136b) 

E=aN J 

Kantrowitz. — Kantrowitz’ theory leads to the third-order 
differential equation 


E=aN+K a J 

Moreland advanced theory. — Moreland’s advanced theory 
leads to the fourth-order differential equation 

E D^yaEE D 3 yo~\~E D^yo'hE Ey o -\-E 1 y o =0 (132a) 


where 


E = E’EiV 1 

E=IM 

E=Et?+a*TNv - 
E=a i N+aTNv 


(132b) 


E^aN 


Moreland elementary theory. — Moreland’s elementary 
theory gives the differential equation 

a Dy o +y o =0 (133) 

Taylor. — Taylor’s theory leads to the third-order differ- 
ential equation 

E IPy o+ E Dhj 0 + E I>Vo+ Eyo=0 (134a) 


where 

E=I »*■ 1 


Ei = 0 
- E=Ea* 
E=K\a 


(134b) 


Temple elementary theory. — The second-order differential ' 
equation for the Temple elementary theory is 

E lEya+E Dy 0 + Egj a =§ (135a) ' 


where 


E £ >3 yo+ E -DV0+ E Eyo = 0 

E=WL 


Ei=IiH~am+kK a 

E=o 

E=aN+K a 


(137a) 


(137b) 


Wylie. — Wylie’s theory leads to the third-brder differential 
equation 


where 


E D*Vo+E &V&E Dy 0 + EVo=0 
E a —Ii?L 

E = IiH- akN+ kK a 

E=a i N+aK cc 

E=(aN+K a ) (l+£) 


(138a) 


(138b) 


stability of motion 


The basic equations of motion having been established 
according to the various theories for a rigid swiveling landing 
gear, attention will be directed next to the meaning of these 
different equations with respect to prediction of the shimmy 
behavior of the landing gear. However, before going into 
this subject in detail it may be useful to discuss briefly what 
sort of information is desired about the motion of a landing 
gear. Basically, the most important question is whether 
or not the motion is stable — that is, does the wheel tend to 
move in a straight line (with decaying shimmy oscillations 
or decaying aperiodic motion) or does the tire tend to move 
laterally out from its rectilinear course (with divergent 
shimmy oscillations or divergent aperiodic motion). To 
answer this question of stability for linear systems, the 
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analytic methods of Routh (ref. 26) or Hurwitz (ref. 27) or 
graphical methods similar to those introduced by Nyquist 
(refs. 28 to 31) are available. A brief discussion of these 
methods is given in appendix C. Any of these methods-will 
provide for most cases a procedure for determining whether 
any particular combination of landing-gear parameters and 
rolling velocity is stable or unstable. 

In general, for complicated problems, rather than inves- 
tigate the stability of a landing gear by these methods for 
all possible conditions it may be more convenient and some- 
times more valuable to draw various types of stability dia- 
grams describing the system in question. For example, for 
case I a typical experimental type of stability diagram is 
shown in figure 7, which presents boundaries between the 
regions of stability and instability as functions of trail and 
rolling velocity for a specific landing-gear model. Another 
useful type of stability diagram for some problems might 
be a plot of boundaries between stable and unstable regions 
as functions of damping moment and rolling velocity. 

In determining these stability boundaries, use is made of 
the well-known fact that the motion of a linear system can 
change from a stable to an unstable condition" only where 
the motion is purely oscillatory — that is, where any variable, 
for example is of the form 


(139) 

or where the motion is purely uniform, of the form 


(140) 

Thus all possible stability boundaries can be obtained by 
directly substituting expressions of the form of equations 
(139) and (140) into the basic differential equations. Tn 
connection with the question of what form of the differential 
equation to use, it is of some importance to note' that the 
final form, where the equation is expressed in terms of one 
variable, is often not the most convenient form to use. For 
example, for case I the purely oscillatory boundaries are 
most advantageously obtained by using equations (115) and 
(117) with the substitutions 


>£= 'Pme 1 ’ 1 * 


yo=yo max e t ^ x+a ^ 


(141) 



Figure 7. — Experimental stability boundary for a tire 29 cm in 
diameter (from ref. 17). 


The advantage of this particular choice arises from the fact 
that it leads to two algebraic equations, one of which does 
not include the damping parameter g. This isolation of the 
parameter g usually" eases' slightly the mathematical labor 
of solving for the purely oscillatory boundaries. 

The equations governing the stability boundaries for case I 
for the general theory and for the systematic approximations 
. are listed in appendix D. 

COMPARISON AND EVALUATION OF THE SUMMARY THEORY AND ITS 
SYSTEMATIC APPROXIMATIONS 

The dual object of the present section is (1) to assess fur- 
ther the value of the summary theory by comparisons 
between the predictions of this theory and the available 
experimental data for case I conditions and (2) to determine, 
by comparison of the relative predictions of the summary 
theory and its systematic approximations, what is tho 
simplest satisfactory systematic approximation to the sum- 
mary theory. Discussion of the previously published theo- 
ries, as applied to case I conditions, is contained in a subse- 
quent section. 

For convenience the following discussion is divided into 
separate considerations of stability-boundary conditions and 
unstable shimmy conditions. 

Stability-boundary conditions. — The present section dis- 
cusses theoretical and experimental stability-boundary con- 
ditions insofar as they are influenced by the tire parameters 
In (where n= 1, 2, . . .), £, K y , c x , N, and r and by hysteresis 
effects. In the major part of this discussion the stability 
boundaries to be considered will be of the type obtained by 
plotting curves of trail against rolling velocity for those 
trail conditions that separate regions of stability and in- 
stability. The general shapes of these stability boundaries 
for case I, according to the summary theory and tho sys- 
tematic approximation theories A to D3, are sketched in 
figure 8 for the special condition of no damping or gyroscopic 
moments (y=r=0). It is seen that the summary theoiy 
and approximations A to C2 each predict that at high speeds 
the motion is stable for large trails and unstable for small 
trails; the horizontal boundary fine is the same for each case, 
and is generally located at a trail roughly equal to the tho 
radius. (This boundary is theoretically completely inde- 
pendent of the spring restoring moment pD,<p and is relatively 
independent of swivel-axis inclination k.) Approximations 
Dl, D2, and D3 fail to predict this boundary. Also, these 
three approximations, together with approximations Cl and 
C2, fail to predict any effect of rolling velocity on the low- 
speed stability boundaries, whereas, according to approxi- 
mation B, for sufficiently small speeds the motion becomes 
stable for all small trails and, according to the higher theories, 
for most of the small-trail region. Also, at low speeds and 
large (usually impractical) trails the higher theories (B and 
above) indicate that the motion becomes unstable at suffi- 
ciently small speeds. The effects of the omitted damper 
and gyroscopic-moment terms would be to reduce the size 
of the regions of instability. 

(a) Effect of higher l n terms: As a first test of the sum- 
mary theory and its systematic approximations there are 
available the experimental data of Yon Schlippe and Dietrich 
(refs. 3, 4, and 5) which were obtained with a small model 
landing gear equipped with a pneumatic tire of 26 cm 
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Velocity, v Velocity, y 

Summary theory and Approximation B 

approximation A 



Velocity, y Velocity, y Velocity, v 


Approximations Cl and C2 Approximations Dl and 03 Approximation 02 

Floras 8. — Qualitative comparison of the stability-boundary predic- 
tions for case I according to the summary and systematic approxi- 
mation theories neglecting damper and gyroscopic moments 
(g=T—0). 

(10 in.) diameter. This model landing gear was tested at 
relatively low speed conditions where the higher Z„ terms 
(/ 2 , 1 - 3 , . . .) are of some importance; consequently, these 
data provide an opportunity for testing the relative and 
absolute validity of the summary theory and the higher 
approximations A to C2 (which differ essentially only by 
then- inclusion or omission of the higher l„ terms.) 

The basic landing gear and tire . constants for the Yon 
Scblippe-Dietrich model, which was tested only in the un- 
tilted condition (k= 0), as taken from references 3 and 5, are 
as follows: 

K—p=g=0 

7* » 0.53 + 0.0025a 2 cm-kg-sec s 
L= 10 cm 
77=640 kg/radian 
77*= 3, 040 cm-kg/radian 
K\— 45 kg/cm 

The quantities l u h, and the higher l x ’s were calculated from 
the previously discussed relations h=li—L, and 

l n =(nL-\-h)h n ~ l /nl. (See eqs. (77) and (19a).) 

The experimental data obtained by Von Schlippe and 
Dietrich for this model are shown in figures 9 and 10, to- 
gether with the corresponding predictions of the summary 
thcoiy and the systematic approximations A to C2. (Also 
shown are the predictions of the theory of Yon Schlippe and 
Dietrich which are discussed in a later section.) Figure 9 
presents stability-boundary plots of trail against velocity, 
and figure 10 presents the frequency at these stability bound- 
aries as a function of velocity. No theoretical curves are 
shown on these figures for approximations Dl, D2, and D3 
since these approximations are too crude to give any de- 
tailed information for this problem; they predict either com- 
pletely stable or completely unstable motion for all positive 
trails (see fig. 8). The equations used to calculate the 
theoretical curves in these two figures are given in appendix 



O 200 400 600 800 1DOO 1200 


Velocity, v , cm/sec 

Figtob 9.- — Comparison of theoretical and experimental predictions of 
the stability boundaries for the Von Schlippe-Dietrich test model of 
references 3 to 5. (Hysteresis effects are neglected in these theoretical 
curves.) 



0 200 , 400 600 800 1,000 1,200 

Velocity, v, cm/sec 


Figtob 10. — Comparison of theoretical and experimental shimmy 
frequencies on the stability boundary for the Von Schlippe-Dietrioh 
test model of references 3 to 5. (Hysteresis effects aro neglected 
in these theoretical curves.) 

D. In these calculations the gyroscopic torque term in- 
volving r has been neglected since r is unknown for these 
data. A rough value of r could perhaps be estimated, but 
such a dubious estimate did not appear necessary because 
the term involving r, according to any reasonable estimate 
of r, would be of no importance in the velocity range of these 
experimental data. 

In comparing the theoretical curves in figures 9 and 10 it 
is observed that approximation A gives a boundary very 
close to that of the summary theory. Approximation B 
does not give as close agreement but it is still fairly good 
and, more important, for most of the trail range the difference 
between approximation B and the summary theory is small 
as compared with the difference between the summary 
theory and the experimental data. As was previously noted, 
approximations Cl and C2 (which axe identical for the 
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present condition of *=0) predict a trad-velocity stability- 
boundary which is independent of velocity so that this 
approximation is an inadequate representation of the 
summary theory at low velocities. However, at high speeds 
approximations Cl and C2 give the same stability boundary 
and frequency as the higher approximations. 

As a further aid in comparing the different systematic 
approximations -with the s ummar y theory, figure 11 presents 
a plot of the linear damping coefficient g required to stabilize 
the motion of the Yon Schlippe-Dietrich model at a medium 
trail of 7 cm as calculated according to the summary theory 
and the various systematic approximations (the equations 
used are presented in appendix D). This figure confirms 
the conclusions drawn from figures 9 and 10 that approxi- 
mation A is a very good representation of the summary 
theory and that approximation B is also a good representa- 
tion of the summary theory. However, more importantly, 
figure 11 demonstrates that approximations Cl and C2 
also give a fairly good representation of the summary theory 
with respect to prediction of the maximum amount of damp- 
ing (i. e., the maximum value of g ) required for stabilizing the 
motion. Approximations Dl, D2, and D3 are seen to give 
inadequate representations of the summary theory. 

The preceding conclusions are, of course, only proven to 
be valid for the specific conditions of the Von Schlippe- 
Dietrich model tests. However, they are believed to be 
valid for most practical rolling conditions. 

In considering the correlation between theory and experi- 
ment for the Von Schlippe-Dietrich test conditions, it is 
noted that the experimental stability boundary in figure 9 
is of the same general shape as that given by the summary 
theory and approximations A and B but that it lies to the 
right of the theoretical curves and thus indicates that the 
experimental system is more stable than the theoretical 
system. Similarly, the experimental frequency-velocity 



Figure 11. — Theoretical calculations of the damping required to 
stabilise the motion of the Von Schlippe-Dietrich model landing gear 
at a trail of 7 cm. (No damping required according to approxima- 
tions Dl and D3.) 


curve in figure 10 falls below the theoretical curves. These 
discrepancies are perhaps a result of the neglect of hysteresis 
damping in the calculation of the theoretical curves of figures 
9 and 10 and are discussed more fully in the next section. 

(b) Effect of hysteresis: In order to investigate whether 
the discrepancies between theory and experiment shown in 
figures 9 and 10 might be explained by a consideration of 
hysteresis effects, some results of calculations involving 
hysteresis effects are shown in figures 12 and 13 for the same 
test conditions as in figures 9 and 10. Figures 12 and 13 
present theoretical calculations of trail-velocity and fre- 
quency-velocity stability boundaries, respectively, both with 
and without consideration of hysteresis effects, together 
with Von Schlippe and Dietrich’s experimental data. Tho 
curves for the Von Schlippe-Dietrich theory were calculated 
by using Von Schlippe and Dietrich’s theory with their 
hysteresis equations (see eqs. 128), whereas the curves for 



Figure 12. — Effects of hysteresis on the stability boundaries for tho 
Von Schlippe-Dietrich test model of references 3 to 6. 



Figure 13. — Effects of hysteresis on the shimmy frequency for tho 
Von Schlippe-Dietrioh test model of references 3 to 6. 
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approximation B were calculated by using the equations 
for hysteresis effects derived in the present report (see appen- 
dix B). The values of and ri a indicated in these figures, 
which were used in these calculations, were estimated from 
static hysteresis-loop data, partly unpublished and partly 
published in reference 24. The values of di and d s used were 
taken from reference 5. It is observed from these figures 
that there is a considerable difference between the calcula- 
tions according to either theory, depending on whether 
hysteresis effects are included or omitted. Moreover, it is 
seen that either theory provides fairly good agreement with 
the experimental data when the hysteresis effects are taken 
into account. These facts suggest that hysteresis effects 
may be an important factor which must be taken into account 
in order to obtain good agreement between theory and experi- 
ment. On the other hand, neither the Von Schlippe-Dietrich 
nor the present report’s consideration of the hysteresis effects 
seems to rest on a completely sound foundation. Thus it 
appeals safe to conclude only that since two different ap- 
proaches to the hysteresis problem indicate that hysteresis 
effects are important, a more rigorous analysis of the 
hysteresis problem would be worthwhile. 

(c) Effect of lii The next test of the summary theory makes 
use of the experimental data of Melzer (ref. 10), who per- 
formed a series of tests with an untilted (k=0) solid rubber 
tire 7 cm (3 in.) in diameter at sufficiently high speeds so 
that his data would be expected to fall considerably to the 
right of the curved-line low-speed stability-boundary curve 
in the first two sketches of figure 8. Eor this velocity range 
the predictions of the summary theory and approximations 
A to C2 are identical; they all predict that the undamped 
(g=T—T\—T a = 0) motion is stable for all values of trail a 
greater than li—N/Ky, that is, the critical value of trail a c at 
which the motion changes from an unstable state to a stable 
one is given by the relation a e =l x . (Even if gyroscopic 
moments are taken into account, this equation is only 
slightly modified throughout a relatively largerange of rolling 
velocity; however, for very large velocities this relation 
breaks down as a result of the gyroscopic moment and the 
motion becomes stable for all positive trails. This phenom- 
enon will be discussed in a subsequent section.) The predic- 
tion a c —li is well confirmed by Melzer’s tests, as is illustrated 
in tho following table of data taken from reference 10, which 
lists the experimental values of h together with the trail 
required for stability for several conditions of vertical 
loading: 


■FV kg 

1 

2 

2.8 

3.0 

Ji=W/AT lt cm 

2.6S 

2.85 

3.22 

3.40 


2.5 

3.0 

3.1 

3.4 



These particular data were taken for the case where no 
spring restoring force acted on the model landing gear (p=0). 
Similar good agreement was obtained for the case where a 
strong linear spring restoring moment was present (p>i£a) . 
For this case, according to the summary theory and approxi- 
mations A to C2, the stability boundary for the positive trail 
condition is the same as for the case of no spring restoring 
force. This prediction is well confirmed by Melzer’s tests, as 
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is illustrated in the following table of data also taken from 
reference 10. 


F,ig . 

1 

2 

2.8 

3.6 

h^NJK x , cm. 

1 96 

2.60 

2.90 

2.91 

' cm — 

2.0 

2-2.25 

• «2.7 ' 

«2.8 


(The difference between the values of Z, in these two tables 
is merely a consequence of changes in tire characteristics 
between the corresponding tests.) 

In order to assess the significance of the preceding compari- 
sons of theory and experiment it should be noted that the 
theoretical relation a 0 =Zi, calculated for g=r=T\=T a = 0 , 
is independent of spring restoring moment p and tire tor- 
sional stiffness K a and, in the velocity range discussed, can 
be shown to be not strongly influenced by gyroscopic or 
hysteresis moments (r, T\, and T a ) or by the higher Z n ’s 
(Zj, Z 3 , . . .). Consequently, this comparison tells practi- 
cally nothing about the correctness of the manner in which 
these important quantities have been inserted into the 
summary theory. On the other hand, the theoretical rela- 
tion a e =l x depends almost entirely on the correctness of the 
kinematic equation of the summary theory which, for 
Melzer’s test conditions, reduces to the kinematic equation 
of approximation C2 (eq. (85)): 

\=—ha 

Thus, the results of the preceding comparison indicate that 
there exists a range of rolling speeds in which the kinematic 
equation of the summary theory, as well as of approximations 
A to C2, is reasonably correct (except possibly for the terms 
involving £, which are as yet not evaluated and are not very 
important). 

It can be said with safety that the range of velocity for 
which the theory gives good agreement with Melzer’s model 
data corresponds to full-scale conditions somewhere inside 
the practical rolling speed range and possibly covering much 
of the practical range. However, the preceding comparison 
definitely does not prove anything about the adequacy of the 
s umm ary theory for small velocities or for the highest 
velocities which may be encountered in practice. 

Further confirmation of the preceding conclusions is pro- 
vided by the experimental data of Scbrode (ref. 17) who per- 
formed tests, similar to those of Melzer, for realistic pneu- 
matic tires as large as 39 cm (15 in.) in diameter, as com- 
pared with the small 7 cm (3 in.) solid rubber tire tested by 
Melzer. Schrode obtained trail-velocity stability-boundary 
plots of the type illustrated in figure 7. These stability- 
boundary plots indicate the same result as Melzer’s data, - 
namely, that there exists a range of velocity in which the 
motion is stable above a certain critical trail a e and unstable 
below it. (This velocity range for the data in figure 7 is 
approximately 60 to 160 km/hr.) It is not possible to check 
quantitatively the theoretical stability-boundary equation 
a 0 =li for Schrode’s* data because Schrode provides no in- 
formation suitable for accurately evaluating Z x . However, 
some qualitative confirmation may be found, since the quan- 
tity Zi always appears to be of the order of magnitude of the 
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tire radius r and for Schrode’s data a 0 is found to be of this 
same order of magnitude (e. g., see fig. 7). Thus Schrode’s 
experimental data appear to confirm the previously drawn 
conclusion that there exists a velocity range in which the 
kinematic equations of the summary theory and approxima- 
tions A to C2 are valid. 

Dietz and Harhng have presented some similar stability- 
boundary curves in reference 16 which also confirm the fore- 
going conclusions. 

(d) Effect of £: Some insight into the effect of the tilt 
parameter £ can be obtained by an examination of the effects 
of swivel-axis inclination k on the stability boundaries ac- 
cording to the predictions of approximation Cl for the con- 
dition where damping, spring-restoring, and gyroscopic 
moments are neglected (g=p=T= 0) in order to isolate the 
effects of inclination. (These assumptions appear to be 
justified for the experimental conditions discussed in this 
section.) Under these assumptions one theoretical stability 
boundary is given by the equation y 


. . %Lh . 

dc—h cos k-\ sm k 

r 


(142) 


Experimental data suitable for testing this relation are avail- 
able in reference 16 for an inclination range of — 20°<A<[20 o 
for one constant-velocity condition. These experimental 
data, some of which had to be slightly extrapolated from the 
data in reference 16, are shown in figure 14 together with 
the predictions of equation (142) for values of £ equal to 0 
and 1. Although reference 16 does not supply the values 
of L, h, and li needed for calculations, the assumed values 
indicated on the figure are probably accurate enough to 
justify the following more or less qualitative observations. 
(The value of l t was chosen so as to make the calculated and 
experimental values agree for the case k= 0.) It is noted 
that the experimental variations and the theoretical variations 
for £=0 are in fairly good agreement and also that these two 
variations are more or less symmetrical with respect to 
positive and negative values of k. On the other hand, the 
theoretical curves for £>0, such as the indicated curve for 
£=1, will all be unsymmetrical. Thus, it appears that £ is 
probably close to zero. In this connection it might be 


-6b 





2b 


o Portly extrapolated experimental data from reference 16. 

} Equation (142); assuming £*5cm, I|=7.5cm, . 

h=2.5cm 


-30 -20 -10 0 19 

Inclination of swivel axis, x, deg 


20 


30 


Figure 14. — Influence of swivel-axis inclination on the stability 
boundary. Tire approximately 12 cm in diameter; jF,= 6.25 kg; 
tt= 19 km/hr. 


noted that Greidanus’ theory, which is the only known theory 
that uses a term corresponding to £, implies a value of £> 1 
(compare eqs. flOO) and (101)). 

(e) Effects of K y and c \ : In order partly to assess the im- 
portance of the tire parameters K y and ci, the special case of 
approximation C2 where r=p=0 and h cos k— a=0 will bo 
considered. While this particular case is of no practical 
importance in itself, its examination permits some insight 
into the effects of the tire parameters K- t and C\. For this 
case the damping required to stabilize the motion is given 
by the inequality 

[a(aN+K« cos /c)-bpJt] gv^+I+gF+gWh^-O 

where p, is given by equation (115c). It is evident from this 
relation that p £ is a stabilizing term if positive and destabiliz- 
ing if negative. Also, according to reference 2, K-, may bo 
as large as F t and, according to reference 24, C\ »0.75. , Thus, 
according to equation (115c), if the .small sin 2 K term is neg- 
lected, p„ may be positive for positive k, whereas if the K y> 
c\, and c Y terms are neglected (as has been done in all pre- 
vious investigations except ref. 2) then p, is always negativo 
for positive k. Therefore, if situations should arise where- 
in p* is important, it is not necessarily safe (o neglect the 
terms involving cx and K y that are used in the determination 
of p,. (See eq. (116c).) 

(f) Effect of cornering power N: As a rough check on the 
variation of the tire cornering power N under dynamic con- 
ditions, there are available experimental frequency data ob- 
tained by Melzer in connection with his previously men- 
tioned tests on an uninclined (k=0) model landing gear 
equipped with a solid rubber tire of 7 cm (3 in.) diameter 
(ref. 10). For the higher velocity conditions of Melzer’s 
tests, the predictions of the summary theory and approxi- 
mations A to Dl lead to the frequency equation 

, • <“*> 

for an uninclined and undamped landing gear, that is, for 
tc=r=( 7 = 0 . (Inclusion of the effect of finite r in this equa- 
tion would not significantly alter this equation for the test 
conditions to be discussed herein.) Some of Molzer’s ex- 
perimental data are compared with the predictions of this 
equation in the following table for the condition p=0. The 


F„ kg 

2.8 

3.6 | 

a/Ii 

0.47' 

0.78 

0.44 

173 

0.88 

f 0p3 

3.8 

KOI 

4.0 

4.8 

ss 

ftrptrimtmtmlt Cp3- --- 

13 

3.6 


b 

B 


experimental data shown represent Molzer’s data for tho 
highest velocity conditions tested. The theoretical and ex- 
perimental values are seen to be in fair agreement. How- 
ever, the experimental values are somewhat smaller than tho 
corresponding theoretical values. This discrepancy is be- 
lieved to be largely due to the fact that these experimental 
tests were not conducted at sufficiently small values of 
shimmy amplitude for the assumptions of a linearized theory 
to be valid. Specifically, all of Melzer's frequency data 
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were obtained for maximum swivel angles of 5° or larger. 
(The data shown in the preceding table correspond to the 
condition d m =5°.) Moreover, Melzer’s data indicate that 
there is a fairly definite decrease in shimmy frequency with 
increasing maximum swivel angle. A sample plot of Mel- 
zor’s data illustrating this effect is given in figure 15. Also 
shown is the theoretical calculation, which is valid only for 
0 m — 0°. If allowance is made for a certain amount of 
experimental error, extrapolation of the experimental data 
to 0 m = 0° could be considered as confirmation of the theory. 
It should be noted, however, that much of the rest of Mel- 
zer’s data, while not necessarily disputing this conclusion, 
do not so clearly support it. Also, plots of the type of 
figure 15 are of limited significance since each test point 
shown corresponds to a different rolling velocity. In view 
of these considerations, the only reasonable conclusion that 
can be reached appears to be that Melzer’s data roughly 
confirms the theoretical frequency and does not conclusively 
dispute its quantitative accuracy. 

Melzer has also conducted frequency tests on the same 
model w r ith an additional strong restoring spring (spring 
stiffness several times the tire torsional stiffness). A com- 
parison of theoretical and experimental frequencies for this 
test is shown in the following table: 


F*t kg 

2.0 

2.8 

3.6 

alh - 

0.77 

0.09 

0.80 

0.69 

0.86 

ftaleulaltd, CPS 

5.2 

6.4 

6.7 

6.6 

6.8 

f *xper{mt*tmlt CpS- -- 

4.9 

6.46 

6.9 

5.8 

5.85 


The much better agreement obtained for this case is ex- 
plained by the predominant influence of the spring restoring 
moment, since for large values of p the model system ap- 
proaches the condition of a simple torsional oscillator with 
moment of inertia 7* and spring constant p, for which con- 
dition the well-known frequency equation is 2n f=-^p/I^. 

In order to assess the significance of the preceding com- 
parisons, the quantities involved in the theoretical equation 
(143) — a, N, JCr, p, and 7^ — will be considered. The quan- 
tities a, p, and 7* are easily measured constants and, for most, 
of Melzer’s data, K a is much smaller than aN; therefore, 
the preceding fair agreement between theory and experiment 
indicates that the quantity N (the tire cornering power), 
which was considered to be a constant in the preceding cal- 
culations, actually does not vary extremely with rolling 
velocity and slummy frequency — at least not for Melzer’s 
test conditions. 


S» 

3 

£ 

ib 


O 10 20 30 40 50 

Moximum swivel angle, B m 

Figure 15. — Influence of shimmy amplitude on the shimmy frequency. 


% 


o Theoretical calculation for 8 m * 0 
a Experimental data from reference 10 


(g) Effect of gyroscopic torque: The next question to bo 
considered is the influence of the gyroscopic torque resulting 
from lateral distortion of the tire. All pertinent experi- 
mental data obtained at very high speeds (e. g., ref. 10; see 
also fig. 7) demonstrate that at sufficiently high speeds the 
previously discussed conclusion that high-speed motion is 
unstable for trails less than h is no longer valid. Instead, 
the experimental data show that at these very high speeds 
instability at any given positive trail ceases above a certain 
critical velocity. The existence of this critical velocity will 
now be shown to result, at least in part, from the gyroscopic 
action which was previously included only in Kantrowitz’ 
theory (ref. 8) but was not specifically mentioned there. 
The simplest systematic approximation that adequately pro- 
vides for this effect is approximation C2. In order to isolate 
the gyroscopic effect, the special condition of no tilt (k=0) 
and no spring-restoring force (p=0) or damper (g=0) will be 
considered. For this condition the equation for the sta- 
bility boundary of approximation C2 (or. Cl) reads 

W+arvX) (a 2 N+aK a +ljTvD =I f v c 2 h ( aN+K a ) (144) 

where the underlined terms are the gyroscopic terms. For 
the computation of the critical velocity v c this equation may 
be simplified still further if it is realized that the quantity 
arli is-small in comparison with the moment of inertia 1$ 
about the swivel axis; hence, for an approximate calculation 
the term artfh can be omitted. . Then solution of equation 
(144) for the critical velocity v c above which the system is 
stable yields the expression ■ 


l(h-a)(aN+K*j 

V hr 


(145) 


which is observed to give an infinite critical velocity for 
zero gyroscopic action (r=0). 

The only available experimental data containing enough 
information on the tire constants that are necessary for 
checking the validity of equation (145) are Melzer’s data 
(ref. 10) and even these data do not provide the required 
gyroscopic moment; therefore, it can only be crudely esti- 
mated as follows: The mass of the tire will be of the order 
of magnitude Wi[2tt ( r—r t )]Trr*, where r is the overall tire 
radius, u the radius of the cross section of the tire torus, 
and w t the average tire density. The moment of inertia 
will be the mass times the radius of gyration r g squared; 
thus, with t!=X according to Kantrowitz, r (eq. (33)) 
becomes 

T^WjrtXr—r^rg 2 
* T= r(r+r 3 ) 


For the usual tire r 4 «0.3r, r 3 is slightly smaller than r 
(say r 3 «0.9r), and r t is probably around 0.8r. Then, to a 
crude approximation, rasOhlu^r 3 . For Melzer’s solid 
rubber tire r= 3.5 cm and w x is probably about 10~ 6 kg- 
set^/cm 4 (specific gravity of 1); thus r«10 _s kg-seeVcm. 
Critical velocities calculated from equation (145) with this 
value of r are compared in figure 16 with some of Melzer’s 
experimental data for one test condition at various values 
of a/h. The calculated and experimental values of critical 
velocity are seen to be of the same order of magnitude. 
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Figure 16.— Comparison of theoretical and experimental variations 
of critical trail with rolling velocity. r= 3.5 cm; fj=3.22 cm; 
^=12.3 kg; 5T a =5.9 kg-cm; ijx=.0;l j7„=0.2. 


Since neglect of the gyroscopic moment gives, theoretically, 
an infinite critical velocity, this agreement indicates that the 
gyroscopic moment is an important factor in producing 
stability at high velocities. Also, the theoretical calculation 
is conservative, that is, the unstable region is overestimated. 
The quantitative agreement between theory and experiment 
is fair but far from excellent. One probable reason for some 
of the disagreement is the relatively crude procedure used 
for estimating the parameter r. Another possible explana- 
tion may he in hysteresis effects, as follows: 

If the differential equation for approximation C2 (see eq. 
(124a)) is modified to take into account hysteresis forces 
and moments in the manner suggested in this report, a 
modified differential equation results (see eq. (B2)), which 
has the stability boundary equation 

+ aW T\v c J r aK a T a v c ) (eW-fgA. + 

liTVe~-\-aNTxv c -\-K a T a v c ) — I+v c \ (aA T +K a ) (146) 

for the same conditions as the corresponding equation (144), 
namely, K=p—g=0. After neglecting the underlined terms, 
which are relatively small at large velocities, this equation 
can be expressed in the simpler and more easily interpreted 
form 

i . aNT\+K a T a (li—a)(aN+Ko) / 117 x 

P ‘ + TTr (147) 

This equation indicates that the effect of finite hysteresis 
(that is, 0; Taf^O) is to reduce the critical velocity be- 
low what it would be for no hysteresis effect (T\=T a =0). 
This result is also indicated in figure 16, where calculated and 
experimental curves are shown for the previously discussed 
high-speed conditions of Melzer’s model tests. (The values 
of T\ and T a needed for these calculations were obtained 
from equations (60), (62), and (143), by using the previously 
mentioned estimated values of ^=0.1 and i? a =0.2, based 
on static hysteresis loops.) In figure 16 the experimental 


data he between the theoretical curves for “hysteresis con- 
sidered” and “hysteresis not considered.” The theoretical 
calculation that includes hysteresis is extremely unconserv- 
ative. Two conclusions can be drawn from these obser- 
vations. First, if the actual hysteresis effect at high speeds 
is only a fraction of the calculated effect, this fact might 
explain the difference between the experimental curve and 
the theoretical curve that does not include hysteresis. Sec- 
ond, as was previously noted, it is evident that the treat- 
ment of hysteresis effects in the present report is inadequate 
and unconservative at high velocities. 

In concluding this discussion of gyroscopic torquo, it 
should be noted that for the case of a rigid landing gear the 
critical design condition (velocity at which shimmy is most 
intense) occurs at low rolling speeds where the gyroscopic 
moment is insignificant. ‘Thus, the inclusion of this gyro- 
scopic moment in the theory is of somewhat academic interest 
(at least for case I) and it probably could be safely omitted 
in practical design calculations. 

Unstable shimmy conditions. — As a further overall check 
of the summary theory and its systematic approximations, 
the experimental data of Kantrowitz (ref. 8) for unsteady 
shimmy conditions are available. 

The significant features of unsteady shimmy motion are 
the divergence and frequency of the oscillation, which are 
simply the real and imaginary parts of the roots of the 
characteristic algebraic equation corresponding to the differ- 
ential equation in question. Kantrowitz has made measure- 
ments of these quantities for a model tiro of 4-inch diameter 
at inclination angles k of 5° and 20° with corresponding 
trails of about 0.08r and 0.31r, respectively (ref. 8). His 
experimental results for k= 5° are presented in figure 17, 
together with corresponding theoretical calculations made 
according to approximation B, which is the simplest system- 
atic approximation to the summary theory which describes, 
at least qualitatively, the shimmy phenomena throughout 
the, complete range of rolling velocity. The theoretical and 
experimental frequencies are seen to be in fairly good agree- 
ment. The theoretical and experimental divergences are in 
fair qualitative agreement, but the experimental variation 
is sometimes considerably below the corresponding theoreti- 
cal one. This quantitative disagreement may be due to 
several factors. First, hysteresis effects are neglected in 
the theoretical calculations. Although use of the hysteresis 
force and moment equations derived in this report would not 
completely explain the disagreement, it is believed that 
these hysteresis equations are not accurate enough, par- 
ticularly at small trails (a=0.08r for the data in fig. 17), 
to justify the conclusion that the disagreement cannot bo 
explained by hysteresis effects. A second partial explanation 
of the disagreement arises from the fact that the theoretical 
calculations may be based on insufficiently accurate values 
of the necessary tire parameters, since Kantrowitz did not 
make direct measurements of all of the most fundamental 
tire parameters, such as h, a, N, and K a . Specifically, he 
measured only the quantity L, a quantity approximately 
equal to aN cos K+K a cos 3 k for 2 values of k, and the path 
frequency and trail a for kinematic shimmy (shimmy with 
velocity approaching zero). The basic tire parameters used 
for calculating the theoretical curves in figure 17 were 
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Figure 17. — Comparison of theoretical and experimental shimmy 
frequency and divergence for Kantrowitz’ experimental data. 
1 k=5°. 

approximately deduced from these quantities as follows: 
The quantity h was obtained from equation (Dl) after 
setting v equal to 0 and substituting Kantrowitz’ experi- 
mental values of L, v u and a for kinematic shimmy. This 
procedure for determining the quantity h is not necessarily 
accurate since equation (Dl) neglects hysteresis effects. 
The tire deflection needed for calculating the trail was 
estimated from figure 8 of reference 2. The trail was com- 
puted from the tire radius, the tire deflection, and the 
inclination. With the aid of this estimated value of trail 
the tire parameters N and K a can be obtained from Kantro- 
witz’ approximate values for aN cos k +K a cos 2 k. In addition, 
most tire tilt effects were neglected; specifically, £ and p, 
were taken equal to zero. While the foregoing procedure 
will probably give roughly correct values for most of the 
fundamental tire constants, it is believed that the limita- 
tions of this procedure and the neglect of hysteresis effects 
are sufficient to prohibit the making of any strong point out 
of the discrepancies between theory and experiment in figure 
17. Thus, to summarize, it appears that Kantrowitz’ data 
furnishes only a rough overall confirmation of the summary 
theory. Although quantitative agreement is poorer than for 
most of the previously discussed experimental data, this 
poorer agreement is not necessarily significant. 

This completes the discussion of case I with respect to 


the summary theory and its systematic approximations. 
Next, attention will be directed to a discussion of case I 
with respect to the predictions of some of the previously 
published theories. 

PREDICTIONS OF SOME OF THE PREVIOUSLY PUBLISHED THEORIES 

Some interesting features of the previously published 
theories in relation to case I are discussed in the following 
paragraphs. Comments on the influence of swivel-axis 
inclination will not be repeated here. 

The theory of Von Schlippe and Dietrich gives predictions 
which are substantially the same as the predictions of the 
summary theory, as can be seen by a comparison of the 
predictions of these two theories in figures 9 and 10 for Von 
Schlippe and Dietrich’s model test conditions. Hysteresis 
effects were neglected in computing both sets of theoretical 
curves. In comparing these two theories it should be noted 
that the only difference in the two sets of theoretical curves 
arises from a slight difference in the expressions used for the 
elastic forces and moments of the tire (see section entitled 
“Forces and Moments on the Wheel” and the comments 
after eq. (128)). The Yon Schlippe-Dietrich theory also 
provides for some tire width effects, but these effects are 
believed to be relatively small for the present test condi- 
tions and were not taken into account in computing the 
theoretical curves in figures 9 and 10. From these figures 
it is seen that the differences between the stability bound- 
aries and frequency .curves for the Yon Schlippe-Dietrich 
theory and the summary theory are usually small in com- 
parison with the differences between the theoretical curves 
and the experimental data. Thus, it seems reasonable to 
conclude that there is no significant difference between the 
main features of the summary theory and the Yon Schlippe- 
Dietrich theory. 

Bourcier de Carbon’s advanced theory provides essentially 
the same predictions as approximation B and will thus 
probably give a reasonable prediction of shimmy behavior 
for the complete velocity range. Similarly, Bourcier de 
Carbon’s elementary theory, corresponding to approximation 
C2, will probably give reasonable predictions for the high 
velocity range. 

Melzer’s theory correctly predicts the existence of the 
large-trail stability boundary given by the equation a e =l u 
but it also predicts the existence of stable motion in the 
small negative trail region between zero trail and a trail 
equal to —e=—K a /N. The latter prediction is in dis- 
agreement with the experimental data of Yon Schlippe and 
Dietrich (ref. 5) who conducted some tests in this trail range 
and found the motion there to be unstable. 

The stability boundary according to Moreland’s advanced 
theory for the case of no damping or spring-restoring forces 
(see eq. (132a)) is given by the equation 



This equation is plotted in figure 18 for zero time constant 
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Figure 18 . — Variation of critical trail with rolling velocity according 
to Moreland’s advanced theory. 


(for which, case Moreland’s theory reduces to the subcase 
of approximation C2 for which e=K a =0) and for several 
finite values of the time-constant parameter r 2 . It is seen 
that if n is large, a large-trail stability boundary no longer 
exists at the trail a 0 =Zi. Since the actual existence of this 
large-trail stability boundary has already been demonstrated, 
it appears likely that t 3 cannot be very large. On the other 
hand, if r 3 is small the introduction of the time-constant 
term produces an almost linear decrease of critical trail 
with increasing velocity until a certain limiting velocity 
(equal to h/ T) is reached; above this velocity all motion is 
stable. Thus, the influence of the time-lag constant term is 
somewhat, like that of the previously discussed gyroscopic 
moment due to tiro distortion or the hysteresis force and 
moment, which may also produce stability at high velocities. 
However, in regard to the general shape of the critical 
trail-velocity curve, the variations predicted by considera- 
tion of the gyroscopic or hysteresis effects (see solid and 
dashed lines in fig. 16) appear more like those of the experi- 
mental data (figs. 7 or 16) than does the nearly linear 
variation predicted from Moreland’s time-lag term for 
small tj. 

Moreland’s elementary theory, Temple’s elementary 
theory, and Maier’s and Taylor’s theories are too crude to 
give any details for case I. 

Kantrowitz’ theory incorrectly predicts instability for all 
positive trails in the absence of damping or gyroscopic 
moments. (This prediction is a consequence of the fact 
that the third coefficient in Kantrowitz’ differential eq. 
(137a) is zero.) 

Wylie’s theory (see eqs. (138)) correctly predicts the 
existence of stability at large trails; however, the predicted 
value of critical trail is given by the equation 

a e ( a c + e)Nl i =Iv t h (149) 

for k= 0. This relation implies that the critical trail is a 
continuously increasing function of velocity, whereas the 
previously discussed experimental data clearly indicate that 
the critical trail rapidly reaches the maximum value h. 
Also it is noted that the last of Wylie’s equations (138b) con- 
tains the factor 1+a/L, which does not appear in any of 
the other theories. This term appears as a consequence of 
the earlier mentioned fact that Wylie’s theory does not 
correctly predict the influence of trail for the yawed ro llin g 


condition. It might be noted that if Wylie’s theoiy wero 
modified to remove this difficulty, as suggested in equation 
(110), then this factor 1+a/L would be replaced by 1; thus 
the modified Wylie theory would be more consistent with 
the other equations. Also, for high velocities the modified 
Wylie theory would more correctly predict a definite critical 
trail according to the relation a c =L. 

PRACTICAL APPLICATION 

Before concluding the discussion of case I, it is perhaps 
pertinent to comment on whether the preceding theoretical 
predictions for this idealized case can bo applied to practical 
landing-gear problems. Moreland has demonstrated that 
neglect of the torsional and lateral elasticity of the landing- 
gear strut can sometimes lead to false predictions of stability 
in an actually unstable system (ref. 12). For example, for 
too large torsional damping the torsional damper unit 
effectively locks the swiveling structure against torsional 
movement with respect to its point of attachment to the 
landing gear strut, so that torsional movement of the 
swiveling structure is possible only as a consequence of the 
always finite torsional stiffness of the strut. In other 
words, for too large damping the damper unit and strut 
combination behaves like a pure torsion spring and is thus 
ineffective for damping purposes. In regard to lateral strut 
stiffness Moreland has pointed out that for a rigid the, 
corresponding to approximation D3, if the strut is pro- 
hibited from lateral motion, as is assumed for case I, the 
motion is stable; on the other hand, if the strut is assumed 
to be of stiffness approaching infinity, the system may bo 
unstable. While this particular criticism of tho practical 
value of case I for infinite stiffness applies only to approxi- 
mation D3, still for finite strut stiffness of the same ordor 
of magnitude as the the lateral stiffness the theoretical 
equations of the other approximate theories may also bo 
un conservative. 

CASE n 

This section is concerned with the idealized landing gear 
shown in figure 6. This landing gear consists of a wheel 
free to swivel about an uninclined ahvays-vertical swivel axis 
that is attached by a horizontal linear spring, of spring 
constant k, to the supporting structure. This configuration, 
case H, is discussed here for two reasons: it illustrates tho 
effect of structural elasticity on wheel shimmy, and it is 
better suited than case I for evaluating approximations 
Dl, D2, and D3 as applied to landing-gear problems in- 
volving structural elasticity. (It may be recalled that theso 
three approximations were of little value in dealing with a 
rigid landing-gear strut (case I); however, for a flexiblo 
strut these approximations may sometimes be of value.) 
In discussing case II, and also later cases, no further men- 
tion will be made of the previously published theories or of 
the question of agreement between theoiy and experiment; 
all discussion will be restricted to the summary theoiy and 
its systematic approximations. 

The discussion of case II proceeds as follows: Tho equa- 
tions of motion for case II are derived according to tho 
summary theory and its systematic approximations. As 
for case I, it is more convenient to rederive these equations 
of motion in a slightly different manner rather than to qpply 
the equations derived earlier for the completely general 
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case. After these derivations ore made the equations for 
the stability boundaries are established. Finally, the 
damping required to prevent shimmy is presented in curves 
as a function of strut stiffness and rolling velocity for a 
specific sample landing-gear configuration according to the 
predictions of approximations C, Dl, D2, and D3. (For 
the present case approximations Cl and C2 are identical 
and are henceforth referred to collectively as approxima- 
tion C.) These curves are utilized to obtain some insight 
into the accuracies of the predictions of appro xim ations 
Dl, D2, and D3 with respect to the more advanced approxi- 
mation C. 

GENERAL DERIVATION 

The derivation of the equation of motion for the summary 
theory proceeds as follows. The details of the landing 
gear considered are illustrated in figure 6. This gear has a 
rigid symmetrical swiveling part having a mass m and a 
moment of inertia about its center of gravity J 0 . The 
nonswiveling part of the landing gear consists of a spring 
of stiffness k with an attached mass m x . The lateral dis- 
placement of the swivel axis is designated as i) a . 

Setting the sum of the lateral spring and inertia forces 
acting on the swiveling part equal to the inertia reaction 
of its center of gravity m D 2 (t} a —c a d) yields the relation 

D * T] a =m D?i\ a —'m<h D<*0 (150) 

Substitution for Xo from the relation 

Xo=f/o— vo—Vo—Va+ad (151) 

(see fig. 6) yields, after rearrangement, 

K\y 0 — D, 2 +m D, 2 +K\+k) V a+ (m<h D t 2 -\-aKk)d= 0 (152) 

Setting the sum of the moments about the center of 
gravity of the swiveling part equal to the inertia reaction 
yields the result 

K c fl—K\\fii—kr} a c 2 —'m l D, 2 r] a c 2 —g D,6—pd—TV D t \ 

=io D?0 (153) 

(see fig. 6) where I 0 represents the moment of inertia of the 
swiveling structure at its center of gravity (J 0 =/#.— racj 2 ). 
Substituting for a and Xo according to equations (23) and 
(151) then yields, after rearrangement, 

(tv D,~ KaV-'Dt+KxcJyo+imiC-i TD?—tv D t +kc 2 —K^c 1 )rj a + 
(Ia k) 2 -\-g D t -\-Tav D ^ p-\-K a -\-aciK\)d=0 (154) 

The third equation for this system for the general case 
is given by the kinematic relation of equation (19a) or 
(19b). When y is omitted, the space derivatives are re- 
placed by time derivatives, and i? 0 is set equal to y a —a6, 
> this relation becomes 

— (1 -\-l l v- l D t +l i v~ i D<‘+ . . .)y 0 +(h— a)6+ Va =0 (155a) 

from equation (19a) or 

— (1 +Lv~ 1 D,)(?‘ , ~ 1d ‘ y 0 +(h—a)d+iia=0 (155b) 

from equation (19b). 

The three equations (152), (154), and (155a) or (155b) 
completely describe the motion of the landing gear according 


to thq summary theory in terms of the three variables 
y 0 , Va, and 0. The corresponding equations for the system- 
atic approximations are obtained as follows. 

SYSTEMATIC APPROXIMATIONS 

Approximation A. — For case II, the three governing equa- 
tions of motion for approximation A are the force and 
moment equations (152) and (154) and the kinematic 
equation: 0 

— (1 + l y v- l D ,+ 4tr 2 D, 2 + y~ 3 D t 3 )y 0 + (li—a)d+va=0 (156) 

Equation (156) is obtained by omitting all l n ’s for n.)> 3 in 
the general kinematic equation (155a). ' 

Approximation B. — The three governing equations of 
motion for approximation B are the force and moment 
equations (152) and (154) and the kinematic relation 

— (1 +l 1 v- 1 D,+L 2 v~ 2 D l 2 )y l) + (4— a)e+va=0 (157) 

which is obtained by omitting the 4 term in equation (156). 

Approximation C. — The kinematic equation for approxi- 
mation C is obtained by omitting the 4 term in equation 
(157). The resulting relation is 

— (l+Zi»“ 1 D f )yo+(4 — &)0+i?a=O (158a) 

or 

v*=(i 9 (158b) 

The force and moment equations (l62) and (154) also apply 
for this approximation. However, a slightly simpler form 
of these equations can be obtained by substituting for y a , 
according to equation (158b), in the terms containing 
K\-r)a. and rrj a in these two equations. With this substitution, 
the force equation becomes 

—Nv~ l D t yo— (t«i D t 2 +rn D t 2 +k)t] a + (me? D, 2 +N)0=0 (159) 

and (using the relation -fir a =elV) the moment equation 
becomes 

^ yo+( m i c -2 Dt 2 +kc 2 )va+ 

[j 0 D t 2 +g D t +~^+ P +(Ci+t)Nj 0=0 (160) 


Equations (158a), (159) (or (152)), and (160)- (or (154)) 
describe the motion for approximation C. 

Approximation Dl. — The equations of motion for approxi- 
mation Dl are obtained by setting Xo equal to 0, or 


and 


y 0 =va—ad 


(161) 


in the force and moment equations (159) and (160). Thus 

— (m x D t 2 +m D 2 -\-Nv~ 1 D t +k)i] a + 

( mca D t 2 -\-aNv~ 1 Dr(-N)d= 0 (162) 

[mjCj, D t 2 —Nv- 1 (c l + e)D t +kc 1 ,] Va+[Ia D t 2 +g D t + 

cbNv~* (c x -\- e)D t~\- p-\- (ci-(-e)fV] 0=0 (163) 

Equations (162) and (163) completely describe the motion for 
approximation Dl. 
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Approximation D2. — The equations for approximation D2 
are obtained as follows. In the force equation (152) set a 
■equal to 0, or 

0 =® _1 DtVo (164) 

which gives the relation 

(mc i ir 1 Dt a +aI&tr 1 Dt+l&)yo— 

* (m t D t 2 +m D?+K x +k)T , a = 0 (165) 

Tor the moment equation (153), set a equal to — \>/li (eq. 
(85)) according to approximation C, and apply equations 
(151) and (164). The result is 

[Idp- 1 D t 3 +gv- l D, 1 +aTDt 1 +TvD t +pv~ 1 D t +a(c 1 +e)Exv~ 1 D t + 
(ci+e)I£x] y 0 +[»hC2 D?—rv D t —(c 1 +€)K x +cJc] i?«= 0 (166) 

Equations (165) and (166) are the basic equations for 
approximation D2. 

Approximation D3. — The basic equation for approxima- 
tion D3 is obtained by first solving the equations (162) and 
(163) of approximation Dl simultaneously to eliminate 
either ij a or Q and then letting N approach « so that all terms 
not multiplied by N vanish. The resulting equation, after 
dividing out the factor N and using the relation c 1 -\-c i =a, 
can be expressed in terms of 0 as 

{[J 0 +mia(a+e)+mci(ci-l-€)] H, 3 + { [m, («+<;)+ 

t m( Cl +e)] v+g}D*+[ P +aH*+e)] Z?rf ®i(a+ e ))0= 0 {167) 

STABILITY BOUNDARIES 

The stability boundaries for case II are obtained in the 
same manner as those for case I. Eor the summary theory 
they are obtained as follows. 


Purely oscillatory boundaries. — The equations for the 
Purely oscillatory motion boundaries are obtained by substi- 
tuting into the differential equations the expressions 

B=6 aax e i ’ t 

Va=Va v ^e i <’ ,+ ^=i laaax e i ’ , (co9 ffi+i sin <r,) - (168) 

yo = 2/Omax el( ’" + ' 1> =i/»«a I e< ’‘( C03 8in 
Substitution of these relations into equations (162), (154), 
and (155b), differentiation and cancellation of e tr ‘, and 
separation of real and imaginary parts into separate equations 
yields the expressions: 

Erom equation (152), 

%\(yo maT cos <r 2 ) + (m 1 v l -\-mv l —K\—k) (va max cos o-0-f- 

(aKx— mc2^)0 nai =O (169) 

and 

^x(l/o MI sin + (my+r/uP—ICx—k) ( Vania sin 0-0=0 (170) 

Erom equation (154), 

CiKxiy 0aax cos o- 2 ) - {rw-KaV-'v) (y 0jxax sin o- 2 )-f- 

(— cos <n)+TVp(na max 8ULO-0 + 

(-IoS+p+Ka+actKje^O (171) 

and 

CiKx(y 0xaz sin <r 2 ) + (jw—Kjr'v) (ij 0max cos <r 2 )-f 
(— mc^+kcr- K\Ci)(va max sin m)— 

TWfyamax 003 ipv+Taw)d max — 0 (172) 


Erom equation (155b), 

“ 008 ^)+P^(Vo wlax sin o- 2 )-Hfi-a)0* ai +(i7 w cos cri)=0 (173) 

and 

-P2»(l/o m « cos sin ^)+ V a maj: sin or=0 (174) 

Equations (169) to (174) can be considered as six linear simultaneous algebraic equations with no constant terms in the 
five variables Vo max cos <r a , Vo aax sin o- 2 , Va^ cos op Va max sin op and d maz . Then for this system of equations to have solutions 
other than zero it is necessary that the determinant of the coefficients of any group of five of these six equations should equal 
zero. The determinant for equations (169), (170), (171), (173), and (174) is 


Ci K\ 

— TW-\-Ela v h 

— niiCiV 2 + c Jc—Cikx 

tPj> 

— To 1 ' 2 + p 4~ -f- ac{K\ 

Kx 

0 

(m 1 +m)i' s — K\— k 

0 

aKx—moxv 1 

0 

Kx 

0 

(mi+m)^—K\—k 

0 

-Pi 

Pi 

1 

0 

l\ — & 

-Pi 

— Pi 

0 

1 

0 


( 176 ) 
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and the determinant for equations (169), (170), (172), (173), and (174) is 


TW — KaV^V 

c x K, 

— tVv 

— mcn^+kci— K\c x 

gv-\-Taw 

K\ 

0 

( m x -\-m)i?—K\—lc 

0 

aKx—maf 

0 

K\ 

0 

(mj -(-m)F 1 — K\— k 

0 

-Pi 

Pi 

1 

0 

l x —a 

-Pi 

-Pi 

0 

1 

0 


(176) 


where p x —p Xa and p 2 =pt a for the summary theory. The corresponding equations for approximation A are obtained by 
setting pi equal to 1— kvi* and p 2 equal to l x v x —hv X: for approximation B by setting p x equal to 1— and p 2 equal to l x v u 
and for approximation C by setting p x equal to 1 and p 3 equal to l x v x . 


The two equations (176) and (176) completely describe 
the conditions governing purely oscillatory modes of motion 
according to the summary theory' and approximations A 
to C (other groupings by fives of the equations merely lead 
to repetition of these two relations). The procedure for 
obtaining the stability boundaries for each of the other 
systematic approximations (Dl, D2, and D3) is similar to 
that just outlined. The resulting stability-boundary equa- 
tions for the other approximations are fisted in appendix E. 

Purely uniform motion. — For purely uniform motion all 
variables will have constant values which may be represented 
as 

0=6max 

Va—Va aaI 

y°~y°max 

Substitution of these relations into equations (152), (154), 
and (155) yields the results 

aK\d max — (.K\-\-k)r)a aax -\-K^y 0max =0 
i.p J rK a J rdC\K))d max J r (cjc— Cii^)5?a CTai + CiK\yo maj — 0 
(l x d)0 nax -\-ria max J/o raaI =0 

For nonzero solutions of these three equations the determi- 
nant of the coefficients of 0 max , Va max , and Vo max must be zero. 
Evaluation of this determinant gives simply 

a+H-^=° (177) 

CHARACTER OF THE MOTION BETWEEN STABILITY BOUNDARIES 

In order to determine the character of the landing-gear 
motion (stable or unstable) between stability boundaries it 
is first convenient to solve the equations of motion for each 
approximation simultaneously to obtain a single linear 
differential equation in one variable for each approximation. 
From these differential equations the stability of the motion 
may be determined by examining the corresponding char- 
acteristic equations by any of the methods discussed in 


appendix C. These characteristic equations for case II, 
according to the various systematic approximations, are 
fisted in appendix F. 

■ EVALUATION OF APPROXIMATIONS Dl. D2, AND DJ 

In the earlier discussion of case I it was not possible to 
present a fair relative evaluation of the three parallel 
approximate theories Dl, D2, and D3 since for case I none 
of these theories provides any realistic information. How- 
ever, for case II such a comparison can be made between the 
predictions of these three approximations and the more 
accurate approximation C, and a specific example will be 
discussed here for a sample landing-gear configuration having 
the relative dimensions and properties A=0.8r, 7i=a==0.5r, 
Ci=c*=0.25r, «=0.3r, m,=0.35m, and r—p= 0. 

The actual calculated behavior of this landing gear in terms 
of damping required for stability as a function of rolling 
velocity according to approximation C is shown in figure 19 
for four values of the ratio of strut stiffness to tire stiffness 
klK\. It is seen from this figure that as the stiffness of the 



Figure 19. — Influence of strut stiffness on damping required for 
stability according to approximation C for a sample landing gear. 
L=0.8r; A=o=0.5r; ci=ci=0.25r; «=0.3r; mi=0.35m; I 0 =mT 3 ; 

r=p= 0. 
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strut is decreased from infinity the damping requirement is 
increased. Also, for large values of strut stiffness the region 
of maximum damping required lies at low speeds, whereas 
for small values of strut stiffness it lies at higher speeds. 

The theoretic predictions of the three theories Dl, D2, and 
D3 for this sample landing gear are compared in figure 20 
with the corresponding predictions of the more accurate 
approximation C (from fig. 19) for three values of strut 
stiffness, k=0.2K\, 1.0 K\, and 5.0 K\. It is seen that for each 
strut stiffness approximations D2 and D3 provide a con- 
siderable overestimate of the damping required for stability- 
On the other hand, approximation Dl gives results in good 
agreement with those of approximation C for the ratios 
kjK\= 0.2 and 1.0, but this approximation greatly under- 
estimates the damping for the large value of strut stiffness 
k/K*= 5.0. 

In view of the comparisons of figure 20, it appears that 
approximations D2 and D3 will not, in general, give reliable 
quantitative estimates of the damping required for stability. 
It appears that approximation Dl may give reasonable 
results for some cases in which the lateral stiffness of the strut 
does not greatly exceed the lateral stiffness of the tire. Since 
this latter conclusion is based on only one set of landing-gear 
parameters, the degree to which it is valid in general would 
require a more extensive investigation for a range of landing- 
gear properties. 

PRACTICAL APPLICATION 

One limitation on the practical application of the preceding 
equations for case II lies in the assumption that the damping 
is simply proportional to the angular swiveling velocity Dfl. 
As was previously mentioned, Moreland has demonstrated 
that this assumption is sometimes unreliable since it implies 


a neglect of the torsional flexibility of the landing-gear strut, 
which in turn can sometimes lead to a false prediction of 
stability for heavily damped systems (see ref. 11 or 12). 
Thus, for systems in which torsional flexibility of the strut is 
important, it will be necessary to replace the damper unit of 
case II by a damper and spring in series, as has been done by 
Moreland (refs. 11 and 12), where the spring represents the 
strut torsional stiffness. This particular case of a series 
damper-spring unit applied to the landing gear of case II is 
not considered separately in this report; it is, however, 
included in the more general case III to be considered next. 

CASE m 

The next type of landing-gear construction to he con- 
sidered is chosen largely to illustrate the application of the 
summary and approximate theories to more complex problems 
than have previously been considered by now making use of 
transfer-function concepts. This landing gear is assumed to 
be of the same general type as that of case II except for tho 
following generalizations. In case II the lateral deflection 
characteristics of the landing-gear strut were represented by 
a single spring and mass combination, or, more precisely, tho 
force exerted on the swiveling part of the landing gear by tho 
strut was set equal to 

^-(i+wi^ii, (178) 

For case III it is assumed that the strut (or, more generally, 
the supporting structure) is a more complex linear system 
than is a spring-mass combination, so that tho strut force- 
deflection relation of equation (178) can bo generalized to the 
new form 

-F,= — Tio(Dt)rj a (179) 



Figure 20. — Comparison of damping required for stability according to approximations C, Dl, D2, and D3 for a sample landing gear. Z/=0.8r; 

A=a=0.5r; Ci = Ca = 0.25r; 6=0.3r; mi=0.35m; Io=mr 3 ; t=p=0. 





THEORIES FOR TIRE MOTION AND WHEEL SHTMMT 


where Tv> is a function of the differential operator D t which 
represents the transfer function correlating F n and rj a - 3 Sim- 
ilarly, the moment on the swiveling part due to the strut, 
which was previously set equal to 

Me= — (p-^-gD t )9 (180) 

for case II, is now generalized to the form 

Me= — T u (D,)0 (181) 

The differential equations for case III are easily obtained 
from the corresponding differential equations for case II by 
replacing lc-\-m t D? by T 10 and p+gD, by T u . For example, 
for approximation Dl equations (162) and (163) are replaced 
by tire new relations 

— [mD t 2 +Nv~ 1 D,-\-Ti 0 (D l )]va+{mc^D, 3 +aNv~ 1 D t +N)d=0 

(182) 

(c 1 + e )N+T u (D,)]6=0 (183) 

which lead to a characteristic equation of an order depending 
on the order of the T’ s. If the T’s are analytically defined 
functions, the calculation of the stability of the motion and 
the stability boundaries proceeds exactly as for case II. 
However, if the T’s are not analytically defined functions 
(for example, if they are determined by experimental tests) 
a slightly different procedure of the following type is needed. 

In order to determine partial information about T 10 , the 
swiveling part of the landing gear can be removed and the 
remaining strut can be subjected to a periodical lateral 
force F n =F Vjnax e ift (either by calculation or by actual vibra- 
tion tests). The resulting lateral-deflection response of 
the structure will have a certain amplitude and phase shift 
which are given by the relation 

F ' (184) 

which is obtained by substitution of the sinusoidal variation 
for -F n into equation (179). The function 1/Ti 0 (w) is a 
complex function of the circular frequency v, the absolute 
value of which represents the amplitude response and the 
argument of which represents the phase shift; it is generally 
called the frequency-response function of the system. 

Similarly, a frequency-response function is defined for the 
response of the landing-gear strut to torsional moment oscil- 
lations by the relation (see eq'. (181)) 

d= ~wi. n M ° ass) 

1 n W 

With the aid of the experimental or calculated functions 
TioO'r) and T n {iv), the stability boundaries for any of the 
theories may be obtained by the usual procedure of sub- 
stituting expressions of the form e irl into the corresponding 

* Moreland, in reference 11, has advanced a similar generalization of the strut lateral-defleo* 
tlon characteristics by means of a concept of "virtual elasticity." However, Moreland's 
generalUation Is less general In that It docs not provide for the existence of strut structural- 
damping forces. 
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differential equation together with the T’s. For example, 
for approximation D3 the basic differential equation is 

{ [Jo+mc, (c,+ t )]D?+ \ mfo+6) r]Z>, 2 + [T„ (iv) + 

a(a-\-e) T vs (iv)]D t J rv(a-\- 1 ) T iri (iv ) } 0=0 (186) 

(obtained by converting eq. (167) to apply to case EH), and 
the stability boundaries for purely oscillatory motion 
(obtained by setting 0 equal to e ut in eq. (186) and separating 
real and imaginary parts) are given by the simultaneous 
equations 

—m(ci+ t)vS— vI[T n (iv) -\-a(a+ e) T r() (iv)]-{- 

r(a+6)R[T 10 (rV)]=0 (187) 

— Tor 3 — rnc x (Ci+ e) v*-\- rR[T n (wO+tt-ffi+e) T M (iv)\ + 

r(a+e)I[Tio(rV)]=0 (188) 

where R and I represent the real and imaginary parts of the 
bracketed functions. Analogous equations are obtained in 
a similar manner for the other approximations. 

In regard to the question as to whether any particular 
motion between stability boundaries is stable or not, case 
IH may present a more difficult problem if the forms of the 
T-fimctions are not knowtf in terms of ratios of polynomials, 
that is, if only the frequency-response variations are known. 
In this event, for example, the usual form of the Routh- 
Hurwitz stability criteria (which is applicable to poly- 
nomial forms only) cannot bo applied and criteria of the 
Nyquist type must be used. A brief discussion of these 
criteria is contained in appendix C. 

. The procedure for applying the summary and systematic 
approximation theories to cases of arbitrary complexity is 
essentially the same as the procedure discussed above for 
case EH, the only important difference being that for the 
general case the equations of motion (eqs. (16) or (19), 
(63), (64), (65), (67), (69), and (71)) are more numerous 
and more complicated. No new concepts need to be 
discussed. 

CONCLUDING REMARKS 

This report has presented a correlation and evaluation of 
the previously published theories of linearized tire motion 
and wheel shimmy and has demonstrated that the major 
merits of all of these theories are contained in a summary 
theory which represents a minor modification of the basic 
theory of Von Schlippe and Dietrich. In cases where there 
are strong differences between the existing theories and 
the summary theory, the previously published theories have, 
in the main, been demonstrated to possess certain deficiencies 
except for Moreland’s advanced ’theory, for which no ade- 
quate evaluation was possible. 

A series of systematic approximations to the summary 
theory has been developed herein for the treatment of 
problems too simple to merit the use of the complete sum- 
mary theory. These systematic approximations have been 
shown to resemble closely the previously published theories 
except that in some details they avoid some of the limitations 
encountered in these theories. 
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Comparison of the existing experimental data with the 
predictions of the summary theory and its systematic 
approximations has indicated a fairly good degree of cor- 
relation between the higher approximations and the existing 
experimental data for the cases investigated. However, 
since the agreement is far from perfect in some respects and 
since most of the limi ted amount of existing experimental 
data was obtained with small models there still remains the 
question as to whether the theory is safely applicable to 
full-scale conditions. In particular, the importance of 
hysteresis damping remains undetermined. 

In regard to the determination of the various tire constants 
required for theoretical shimmy calculations, it is noted 
that the existing pertinent experimental data, mostly con- 
tained in references 21, 24, and 32 to 37, are extremely limited 
and apply mostly to small, obsolete, or foreign tires. Fur- 
thermore, although various attempts have been made to 
correlate and to reconcile theoretically the experimental 
data (e. g., ref. 2 or 38), there still apparently does not 
exist any fully reliable theoretical means for predicting all 


of the needed elastic characteristics of tires. In view of 
these considerations, a need exists for additional experi- 
mental data on modem tire characteristics and also for a 
more adequate evaluation of the existing data to determine 
whether these data can be applied by scale laws to predict 
the characteristics of any tire with tolerable accuracy. 

In regard to the adequacy of a linearized theory of tire 
motion, it can be stated only that there is as yet no strong 
indication that a nonlinear theory is required for prediction 
of the stability boundaries. If, however, a knowledge of 
the large-angle (nonlinear) behavior is required, a theoretical 
system for dealing with this problem could be developed on 
the basis of assumptions of the type advanced by Kelley, 
Rotta, and Temple (see refs. 18, 19, and 21 (p. 36), respec- 
tively). 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 18, 1956. 
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APPENDIX A 

CALCULATION OF EQUIVALENT VISCOUS DAMPERS 


GENERAL CASE 

In the derivations and equations in this report, only 
linear damping terms were introduced so that the resulting 
equations would remain linearized. However, the damping 
moments caused by friction, hysteresis losses, and the ordi- 
nary shimmy dampers are no n! in par and therefore it is nec- 
essary to replace these nonlinear moments with equivalent 
linear viscous moments. The equivalent viscous moment is 
usually determined by assuming that linear and nonlinear 
damping moments are equivalent if they dissipate the same 
amount of energy during each cycle of shimmy oscillation.' 

For a linear damper of moment gD t p the energy dissipated 
per cycle of sinusoidal oscillation is (for 4'='Pmax sin vt) 

i?= f gD t -pd-p 

and, since D t \p—\p ma x v cos vt and dip=ip max v cos vt dt, 


E=gvp m 


■r 


cos 2 vt d(vt) 


which gives 


E=irgvip m 


Therefore the linear damping constant g is related to the 
energy loss per cycle by the relation 


E 


irvtj/ m 


(Al) 


By using this relation an effective value of linear damping 
constant can be calculated for any nonlinear damper if the 
energy dissipation per cycle is known. 

VELOCITY-DEPENDENT DAMPING 

By using equation (Al) Rotta (ref. 2) has calculated the 
effective damping constant for damping moments of the type 

(A2) 

The first term represents friction damping and the second 
represents fluid damping. The exponent will probably 
always be between 1 and 2. 


Rotta’s calculations proceed essentially as follows: First, 
the total energy dissipated is calculated from the relation 


E- 


Mj, d^ 


where 


-I 

'P='Pmax Sin vt 


(A3) 

(A4) 


After combining equations (Al) to (A4) and integrating, the 
following equation for g is obtained: 


g= 


4M 0 


KWpm 


-(Zn&tt*" 1 1 V n ' 1 - f 1 cos" 1+1 vt d(vt) (A5) 
irj o 

It is seen that the damping constant depends on the ampli- 
tude p max ] for tti>l it is large at small angles due to the first 
term and is large at large angles due to the second term. 
The minimum value of the equivalent damping constant 

^obtained by setting equal to 0^ occurs at the angle 


'Pa= 


M 0 




l/n, 


(A6) 


v"‘(ni— 1)J cos" 1+1 vt d(vt) 

and the corresponding minimum damping constant becomes 

(A7) 


”i-i 

g=M 0 


where 


K=- 


4% 


[J Q W2 cos n ‘ +1 vt d(vt)J ln ' (A8) 


tt(% — 1) 

For the special case of velocity-squared damping (n= 2) 
equation (A7) reduces to the relation 


g — q? 

FRICTION DAMPING 


(A9) 


The equivalent damping constant for friction damping 
(constant moment M 0 ), obtained by setting g n equal to 0 
in equation (A5), is 

g= (Aio) 

TTVpmax 
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APPENDIX B 

DIFFERENTIAL EQUATIONS FOR CASE I WITH INCLUSION OF HYSTERESIS EFFECTS 

If the differential equations of approximations B and C2 for case I (see eqs. (122) and (124), respectively) are modi- 
fied to take into account the hysteresis force and moment expressions derived in this report, the following differential 
equations are obtained: 

For approximation B, 

Eo Dh/o+Si D 3/ y 0 -jrE 2 Dy o -\-Eiy o =0 (Bla) 

where 


E 0 =I,t*l 2 

S)=I^Z 1 + (arc 2 cos ^+go)L+ (a 2 K\T x +K a T a cos 2 k)1 2 v 

£’ 2 =Ir 2 +(a.rr 2 cos K+g D )l l +{.a^+K a cos 2 k+p+p,)Z-i+ (a i N'Tx+a'K a T a cos k)v 
E?, — a?N-\- a 'K a cos k + pZi -)- pj-i -\-gv-\-aliTi? cos 2 K-\-(aaNT\ cos n.-\-K a T a cos 2 k)v 
Ei—craN cos K-\-K a cos 2 K+p+p I -)-'Wt 
For approximation C2, 


where 


Eo TEy^-^E^i Z? 2 y 0 -{--E> Dyo-\-Eqyo=0 


E 0 =I i ,v i l 1 


E l =I 4 # l +(aTV 2 cos K+gv)l 1 +(a 1 NT\+aR a T a cos k)v 

Ei=a 2 N+aK a cos K+pk+pJi+gv+liTt? cos 2 K+(aNT x cos ic+K a T a cos 2 k)v 
Ez—aA T cos K+K a cos 2 K+P+p«+«rt 


(Bib) 


(B2a) 


(B2b) 


J 
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APPENDIX C 

STABILITY CRITERIA 


In this appendix a brief review is presented of most of the 
existing methods for examining the stability of motion of 
systems whose motions are governed by linear differential 
equations of the type 

f(D,)d = 0 (Cl) 

The solution of this type of equation consists of terms of 
the form 

6=e pt (C2) 

whence 

D t n 0=p*e v, =p n e (C3) 

Substitution of equation (C3) into equation (Cl) yields the 
algebraic equation 

/(?)= 0 (C4) 

for the p’s. Equation (C4) is called the characteristic 
equation of the differential equation (Cl). 

The type of motion for the linear system is determined by 
the character of the complex roots of the characteristic 
equation. Most important, the motion is entirely stable if 
and only if the characteristic equation possesses no roots 
having positive real parts. Several procedures are available 
for determining whether a particular characteristic equation 
has such roots with positive real parts. 

One procedure which is useful in cases where the charac- 
teristic equation (C4) can be written in the polynomial form 

a v p n +aip n ~ l + . . . +o„= 0 (C5) 

is the well-known Routh-Hurwitz criterion which makes use 
either of the Routh test functions (ref. 26) or of the equiva- 
lent Hurwitz determinants. In Hurwitz’ form the require- 
ment for stability (or no roots with positive real parts) is 
that Oo and all of the n determinants 

Z)i=d\ 



d! 

< h 

at . 

. . (hv-i 


Ot) 

(h 

ai . 

• • Chu-2 

D v = 

0 

«i 

03 • 

• • 0j o _3 





. . o„ 


for u= 1, 2, ... n must be greater than zero, or, in the alter- 
nate form of Cremer (ref. 39), all of the a n ’s and either the 
even or the odd Hurwitz determinants must be positive. 

This criterion is particularly suited for examining the sta- 
bility of linear systems with polynomial-type differential 
equations of low order. However, for high-order polynomial- 
typo differential equations this criterion may not be the 
480164— B8 18 


easiest to use and for nonpolynomial-type equations the 
criterion is not directly applicable. For such cases use can 
be made of the graphical-type criterion originated by 
Nyquist (refs. 28 to 31). Some discussion of criteria of 
this type is contained in most books dealing with servo- 
mechanisms or feedback amplifiers (for example, refs. 29 and 
31). These references provide the necessary theoretical in- 
formation for applying these criteria and the theory will not 
be repeated here; however, it may be useful to set down 
here, together with an example, one mechanical procedure 
for applying this criterion and a few pertinent comments. 

Consider a differential equation with the characteristic 
equation (B4) for a case where the function j(p) cannot 
necessarily be easily expressed in a simple polynomial form 
which can be handled by the usual Routh-Hurwitz criterion. 
(This may be the case, for example, where part of the func- 
tion j(p) is evaluated from experimental frequency-response 
data.) 

The function f(p) is assumed to be a single-valued func- 
tion of p which is real when p is real. It is also assumed 
that the function f(p) has no poles in the region of the com- 
plex p-plane where the real part of p is greater than zero. 
When the equations of motion are set up in the manner fol- 
lowed in this report, the condition of no poles in this region 
is usually satisfied for actual landing gears since this condi- 
tion implies only that the landing-gear strut, as represented 
by equations or experimental curves, possesses some damp- 
ing or is at least not inherently an unstable structure. For 
example, for the equations of case HI, only the poles of the 
functions T w (p) and T n (p) could lead to such poles. How- 
ever, if, for example, T w (p) had such a pole, equation (161) 
would indicate the possibility of a steady or divergent oscil- 
lating force F, corresponding to this pole, even if the lower 
end of the landing-gear strut were held fixed (%=0); this 
obviously cannot occur in actuality. 

In order to decide whether the motion of the system de- 
scribed by equation (C4) is stable, the following procedure 
may be followed: 

(1) Determine the variation oij(p) in equation (C4) for 
the case of pure sinusoidal oscillations, that is, for p=iv in 
the range 0<y< <= . 

(2) Plot the real part of f(iv) against the imaginary part 
for the complete range 0<y< <» . This will give a curve 
such as is illustrated in figure 21 for a sample case. As v 
varies from 0 to a> this curve will move about the origin 
through a net angle of Jic radians, this angle being con- 
sidered positive in a counterclockwise sense. ( J=1 for the 
case illustrated in fig. 21.) 

(3) Determine the asymptotic behavior of the character- 
istic function j{p) of f(iv) for p-> <» or v~> <*> ; at this limit 
the function will behave osf(p)ccp i or/(jy)ocr*, whence j can 
be determined. 
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(4) Under the preceding restrictions f(p) being a single- 
valued function of p, real when p is real, and having no 
poles in the half-plane E(p)> 0, it can be shown that the 
motion corresponding to the differential equation is stable 
if and only if 2 J=j. (In the sample case of fig. 21, where 
J= 1 and j=2, the motion is therefore stable.) 

In conclusion it might be noted that, although either the 
preceding Nyquist type stability criterion or the Routh- 
Hurwitz criterion can usually be applied to most of the ap- 
proximate equations discussed in this report, they cannot be 
directly applied to some equations of transcendental form 
such as those of the summary theory, since such equations 
may correspond to infinite-order linear differential equations 
(for example, see eq. (120)); hence, an infinite number of 
Hurwitz determinants would have to be evaluated or the 
Nyquist type plot would circle the origin an infinite number 
of times. 



Figure 21. — Variation of the funotion / (t*>) for tho linear differential 
equation having the characteristic equation 

P 4 +6/+23p 3 +50j>+50 n 

JVP> pH-2p+2 “ U 


APPENDIX D 

STABILITY BOUNDARIES FOR CASE I 

The following equations describe the conditions under which purely oscillatory motion is possible for case I for the sum- 
mary theory and the systematic approximations. 

For the summary theory and approximations A to C2, 


and 


^2 (g 2 A\+-S~ (t cos 2 K-bp+pt) (jh^+y^+Kfl-fiA+cxFg sin k)Pi — vipJK. a cos *] (<rlj cos k — a) 

Ii v i\vi 1J rPi)~ TV iT?ii< J h cos K—a) cos k 


(<rli cos k cl) [pi(aK\-\-oJ? t sin «)+>q pi(,K a ~nf) cos *] 

wW+pit 3 ) 


■aro cos k 


where, for the summary theory, 


2h=Pi„=cos vji — Lvi sin vji 


for approximation A, 


Pz—Psm — sin vJi-\-Lvi cos h 


Pi=l — W 


for approximation B, 


Pi=kv i—kfi 


2>i=l— kn 2 


(Dl) 

(D21 


and for approximations Cl and C2, 


Pi=hvi 


Pi=l 

P2=hn 


For approximations Dl and D3 purely oscillatory motion does not exist. 
For approximation D2, 


v>= 


a?K\-)raeK\ cos K+p+pi+rr 2 cos 2 k 


j=v '■ 


cos K+eK\ cos 2 K+ciF t sin k cos k) 


a?K\-{-(uK\ cos ic+p+pt+n? cos 2 * 


-ar cos 


*] 


(D3) 

(D4) 


The stability boundaries for uniform motion are obtained by setting the coefficient of the y 0 terms in the various differ- 
ential equations equal to zero. For example, for the summary theory and approximations A and B the equation 

oaN cos k+Sa cos 5 »c-(-p-}-pt+'U*=0 


describes this stability boundary. 
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APPENDIX E 


STABILITY BOUNDARIES FOR CASE H 


The equations governing purely oscillatory modes of oscillation for case H are as follows for approximations Dl, D2, and D3 : 
For approximation Dl, 


[70+^01(0!+ e)+mia(g+e)]r a —[p+afr(g+e)] 

[(mi+mV— fe] { — [7o(m 0 +m)+mm 1 c./]v 4 +[7 0 7+7n7c 2 :i +m 1 A 7 '(a+e)+m7/(c 1 +e)+p(m I +m)]^— l'[p+77(a+e)] } ' ' 


For approximation D2, 


AT -1 { [m(c l +e)+m l (a+€)]v i ~k{a J ce)}[(m l (L J rmc l )i?+N—ak] 
9 [(m 1 +my-kf+N 1 v-V 


(E2) 


^__[Io(.mi+m)+mmiCi ! ]v i —[Io(K\+k)+p(m 1 +m)+mcj > k+m 1 K},a(a+e)+mK\c 1 (c 1 + e)V+ [p(K- K +k)+kK^a(a+ e)l N 


For approximation D3, 


—v{ {rami+TCimy— [ark-^miKxja^r t)+m.Kx(ci+ e)] v LJ rkK\{a,+ e) } 
r 2 [(m 1 +m)^-(Zx+7)] 

j p_ p-\-ak{g-\- e) 

I 0 +m L a(a+ e)+mCi(c!+ e) 


(E3b) 

(E4a) 


g=v[k(a+ e)v -2 — m(c!+ e)—m l (a+ <e)] 


(E4b) 


APPENDIX F 

CHARACTERISTIC EQUATIONS FOR CASE H 

The characteristic equations for case II are as follows for the summary and systematic approximation theories: 
For the summary theory and approximations A to C2, 


(tv — K j}~ l )p + CiK}, 

-m 

where, for the general theory, 

for approximation A, 

for approximation B, 

and for approximation Cl or C2, 


— (wq+m)# 5 — K\— k 
'miC i ‘p % —-rvp-\-ko 1 —K\Ci 


mc a p t +aK\ 

7 0 p*+ (g -\-Tav )]) + p -\-K a a.CiK\ 


=0 


1 

/(p)=(l+AtrW , ' lp 


l\ Or 


f(j?)=l+hv~ 1 p+kv~ i p i + Uv~ 3 y 
fip ) = 1 ■ +bP~ 1 P+kv~*^ 


/(p)=l+Z 1 p" 1 p 

For approximation Dl the characteristic equation is 


[7 0 («i 1 +m)-fmm 1 c 2 2 ]p 4 +{^(mi+m)+lVff 1 [Jo+mia.(a+e)+mc 1 (c 1 +e)]}y 3 +[(7o+?/ic^)7+p(m 1 -f m)+mA 7 (ci+e)+ 
m 1 2V r (a+ «) +gNv~ 1 }'p % -\- { [a.(a+€)A 7 c _1 4-ff]^+pA r ti _1 }p + [p+ (a+e)IV]&=0 

For approximation D2, 


[7o(m l +m)+wmiC 3 2 ]» J y 6 +[rmia+rmc 1 +£/(m 1 +m)t) _l ]p 4 + [Io(K\ + &)t> 1 +p(m 1 +m)r 1 -\-Tv{m l +m)-\-kmcfo~ l ymia(a- s r6)K\v~ 1 -\- 
wCi(c+e)J7xO _1 ]p'' 5 +[Ta7+^(I7x+7)t) _1 +m(ci+€)A'x+m 1 (a+6)I7x]y 2 +t'rrI'+p(Ax+/:)p _1 +a/:77x(a.+e)p _1 ]^+7Arx(a+e)=0 


For approximation D3 see equation (167). 
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